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ABSTRACT 

As  part  of  a joint  technical  effort  Involving  North  Carolina 
State  University  and  United  Technologies  Research  Center,  an  experi- 
ment was  conducted  to  measure  the  response  of  an  Isolated  turbomachine 
rotor  to  a distortion  in  inlet  axial  velocity.  A once-per-revolution 
sinusoidal  variation  in  axial  velocity  with  an  amplitude  of  approxi- 
mately twenty  percent  of  the  average  axial  velocity  was  generated  by 
an  upstream  screen.  The  response  of  the  rotor  was  studied  using 
pressure  transducers  and  skin  friction  gages  mounted  on  one  of  the 
rotor  blades  and  a velocity  probe  at  the  rotor  exit  plane  as  well  as 
with  standard  stationary  frame  pneumatic  Instrumentation. 

The  rotor  was  operated  in  undistorted  flow  to  establish  the 
quasi-steady  Behavior  of  the  compression  system.  When  the  air  inlet 
angle  was  reduced  past  a certain  limit,  the  rotor  began  to  experience 
rotating  stall.  When  the  rotor  was  operated  in  distorted  flow,  the 
pressures  on  the  surface  of  the  instrumented  blade  were  observed  to 
vary  as  a function  of  the  instantaneous  inlet  angle  * These  variations 
were  greatest  at  the  leading  edge  of  the  airfoil  and  became  smaller 
toward  the  trailing  edge.  This  concentration  of  activity  in  the 
leading  edge  region  is  more  pronounced  than  has  been  observed  for 
isolated  airfoils.  As  the  Instrumented  blade  traversed  the  distortion, 
it  was  observed  to  operate  transiently  at  inlet  angles  below  the 
quasi-steady  stall  point  in  an  apparently  unstalled  condition. 

At  these  lower  inlet  angles,  a separation  bubble  was  found  to  exist 
on  the  leading  edge  of  the  blade. 
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Velocity  measurements  made  at  the  rotor  exit  plane  in  undistorted 
flow  at  moderate  incidence  angles  indicated  that  the  flow  was  essen- 
tially two-dimensional  with  no  radial  flow  being  measured  except  in 
the  blade  wake  regions.  This  was  not  the  case  for  similar  incidence 
angles  in  distorted  flow  in  which  case  significant  radial  flow 
components  existed.  The  circumferential  exit  angle  in  distorted  flow 
was  found  to  lag  the  steady  state  exit  angle.  However,  the  variation 
of  exit  angle  as  a function  of  inlet  angle  in  distorted  flow  did  not 
have  the  same  character  as  for  undistorted  flow.  The  observed 
differences  are  attributed  partly  to  the  three-dimensionality  of  the 
distorted  flow  field  and  partly  to  some  measurement  errors  associated 
with  the  geometry  of  the  probe. 

The  peak  pressure  rise  generated  by  the  rotor  was  lower  for 
distorted  flow  than  for  uniform  inlet  flow  as  should  be  expected. 

However,  at  high  blade  loading  the  pressure  rise  was  less  than  the 
peak  value  for  both  undlstorted  and  distorted  flow,  but  was  greater 
for  distorted  flow  than  for  undistorted  flow.  It  is  believed  that 
the  ability  of  the  blade  to  operate  transiently  beyond  the  steady 
state  stall  point  is  responsible  for  this  behavior. 
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INTRODUCTION 


An  aircraft  propulsion  system  must  be  capable  of  operating 
smoothly  and  efficiently  in  all  parts  of  the  aircrafts's  flight 
envelope.  This  means  that  at  some  point  every  engine  will  be  forced 
to  operate  with  non-uniform  inlet  flow.  For  example,  military 
aircraft  which  are  sometimes  required  to  execute  violent  maneuvers 
may  experience  inlet  distortion  generated  by  the  aircraft  attitude. 
Weapon  firing  is  another  source  of  distortion  peculiar  to  military 
aircraft.  One  source  of  distortion  which  is  common  to  all  aircraft 
is  ground  operation  or  takeoff  in  a crosswind  or  in  gusts.  The 
result  of  operation  with  non-uniform  inlet  flow  may  be  a simple  loss 
of  performance  or,  in  extreme  cases,  the  triggering  of  compressor 
stall.  Since  any  loss  of  thrust  under  the  previously  described 
circumstances  may  have  serious  consequences,  the  engine  designer  is 
faced  with  the  task  of  creating  a compression  system  which  is  insen- 
sitive to  inlet  distortions.  That  is,  the  losa  in  pressure  rise 
associated  with  distortion  must  be  small  and  the  tendency  of  the 
compressor  to  stall  under  the  influence  of  distortion  must  be  mini- 
mized. Furthermore,  the  operating  conditions  at  which  some  type  of 
compressor  stall  might  be  expected  to  occur  must  be  known  ao  that 
they  may  be  avoided.  Finally,  in  the  event  that  the  compressor 
does  become  stalled,  it  must  be  possible  to  recover  from  stall  with  a 
minimum  of  effort  and  as  quickly  as  possible. 
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In  order  to  better  understand  the  inlet  distortion  problem,  a 
more  detailed  discussion  of  compressor  stall  is  in  order.  Compressor 
instabilities  appear  in  two  distinct  forms  known  as  surge  and  rotating 
stall  (Ref.  1).  These  two  pheomena  are  described  in  detail  in  Ref. 

2.  The  description  of  surge  and  rotating  stall  as  presented  in  the 
subsequent  paragraphs  is  essentially  a summary  of  that  presented  in 
Ref.  2.  Surge  is  fundamentally  a timewise  excursion  in  the  mass  flow 
through  the  compressor  simultaneously  involving  the  entire  compressor. 
On  the  other  hand,  rotating  stall  consists  of  one  or  more  regions  of 
severely  retarded  (or  even  reversed)  flow  which  rotate  around  the 
annulus.  These  two  types  of  behavior  are  illustrated  in  Fig.  1.  It 
is  important  to  the  designer  to  be  able  to  predict  which  mode  of 
instability  will  occur  because  the  nature  and  hence  th£  consequencs 
of  these  two  modes  of  instability  are  quite  different. 

It  has  been  stated  that  surge  is  a timewise  excursion  in  the 
mean  mass  flow  through  a compressor.  As  such,  it  is  a relatively  low 
frequency  phenomenon  (frequencies  of  one  or  two  cycles  per  second  are 
typical)  which  involves  large  excursions  in  the  annulus-averaged  mass 
flow.  In  fact,  the  conditions  may  be  such  that  rotating  stall  occurs 
during  part  of  each  surge  cycle.  Surge  can  lead  to  excessive  inlet 
overpressures  resulting  in  damage  to  engine  components.  Since  the 
excursions  in  mass  flow  during  surge  are  large,  it  is  very  possible 
that  the  compressor  will  operate  in  conditions  normally  associated 
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ROTATING  STALL  SURGE 

CIRCUMFERENTIALLY  NONUNIFORM  FLOW  AXIALLY  OSCILLATING  FLOW 


INLET  HOT  WIRE  OUTPUT 


FIG.  1 POSSIBLE  MODES  OF  INSTABILITY  ON  STALL  LINE 


i 


4 


with  un8t al led  flow  over  part  of  the  surge  cycle.  Recovery  from 
surge  might  be  possible  by  changing  the  throttle  setting  or  opening 
internal  bleeds. 

In  contrast,  the  unsteadiness  associated  with  rotating  stall  is 
more  localized.  Once  a rotating  stall  system  becomes  stable,  the 
annulus-averaged  mass  flow  is  constant  with  time.  One  or  more 
regions  of  severely  retarded  flow  develop  and  rotate  around  the 
annulus.  They  may  extend  over  the  entire  span  of  the  blading  or 
cover  only  part  of  the  blades  (for  example,  only  the  hub  region  might 
be  involved).  The  circumferential  extent  also  varies  from  a few 
blade  gaps  to  perhaps  half  of  the  annulus.  The  frequency  of  the 
rotating  stall  phenomenon  is  typically  an  order  of  magnitude  greater 
than  that  for  surge.  A rotational  speed  of  approximately  half  the 
rotor  speed  is  common.  The  variations  in  instantaneous  inlet  angle 
(and  hence  blade  loading)  as  the  stall  cell  or  cells  propagate  around 
the  rotor  can  be  quite  severe.  Prolonged  operation  in  rotating  stall 
may  thus  lead  to  fatigue  failure  in  the  compressor  blading.  Recovery 
from  rotating  stall  can  be  much  more  difficult  than  recovery  from 
surge  because  the  compressor  does  not  operate  in  conditions  normally 
associated  with  unstalled  flow  at  any  time  during  the  occurrence  of 
rotating  stall. 

As  a background  for  the  present  experimental  effort,  the  perfor- 
mance of  a rotor  in  distorted,  but  unatalled,  flow  will  first  be 
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considered.  In  particular,  a circumferential  distortion  in  the  inlet 
axial  velocity  will  be  discussed,  because  this  type  of  distortion  has 
been  found  to  have  more  pronounced  effects  on  rotor  performance. 

This  topic  has  been  extensively  investigated  both  experimentally  and 
analytically  by  a number  of  authors.  A good  sampling  of  various 
approaches  that  have  been  used  may  be  found  in  Ref.  3.  References  4 
and  5 provide  a further  bibliography  of  background  material.  The 
analytical  investigation  of  the  performance  of  a rotor  undergoing 
distortion  must  account  for  the  fact  that,  although  the  distortion 
field  is  steady,  the  blades  themselves  are  subjected  to  an  unsteady 
flow  as  they  traverse  this  distortion.  Their  behavior  may  therefore 
depart  considerably  from  their  undistorted  performance  characteristics, 
especially  when  the  blades  are  heavily  loaded.  As  an  example  of  this 
phenomenon,  the  unsteady  total  pressure  loss  of  a compressor  operating 
in  distorted  flow  is  shown  in  Fig.  2 (taken  from  Ref.  6).  This  curve 
also  points  out  another  important  aspect  of  the  unsteady  flow  field 
in  that  the  rotor  is  seen  to  operate  transiently  in  an  apparently 
unstalled  condition  at  inlet  angles  which  would  normally  result  in 
stall  for  undistorted  inlet  flow. 

Many  methods  have  been  used  to  model  the  blade  rows  of  a turbo- 
machine and  to  include  unsteady  effects.  Adamczyk  (Ref.  7)  modelled 
the  blade  row  as  an  actuator  disk  for  which  the  exit  flow  angle  and  the 
total  pressure  loss  across  the  row  were  given  as  functions  of  the  inlet 
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angle  for  steady  flow.  The  unsteady  effects  were  modelled  as  simple 
time  lags  applied  to  the  steady  characteristics.  Colpin  (Ref.  8) 
extends  this  type  of  analysis  to  compressible  flow.  Both  Adamczyk 
and  Colpin  assumed  that  the  flow  through  the  blading  was  two-dimen- 
sional and  that  the  blading  could  therefore  be  treated  as  a cascade. 

Recent  work  (Ref.  9)  has  applied  the  actuator  disk  approach  including 
three-dimensional  effects.  This  analysis  is  capable  of  explaining 
trends  observed  in  experimental  data  which  could  not  be  explained  by 
two-dimensional  theories.  Three-dimensional  effects  in  swirling 
flows  were  the  subject  of  an  analytical  and  experimental  investigation 
carried  out  by  Greitzer  and  Strand  in  Ref.  10.  This  investigation 
also  led  to  the  conclusion  that  certain  types  of  flows  are  inherently 
three-dimensional  and  can  not  be  handled  using  two-dimensional 

1 

theories.  In  a different  approach,  Mazzawy  (Ref.  6)  modelled  the 
compressor  undergoing  distortion  as  several  compressors  in  parallel 
with  different  inlet  conditions  but  a common  outlet  condition.  These 
few  papers  are  fairly  representative  of  the  analytical  work  being 
done  in  this  field.  To  check  the  validity  of  the  analyses  and,  in 
some  cases,  to  provide  input  data,  it  is  necessary  to  have  a body  of 
experimental  data  representative  of  the  conditions  being  modelled. 

Some  authors  have  conducted  purely  experimental  studies  of  the 
response  of  turbomachine  blading  to  inlet  distortion.  Bruce  and 
Henderson  (Ref.  11)  have  used  a force  balance  to  measure  the  unsteady 
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lift  and  moment  on  a blade  in  an  isolated  rotor  for  sinusoidal 
distortions  in  inlet  axial  velocity.  The  unsteady  pressures  on  the 
surface  of  a blade  were  measured  by  Peacock  (Ref.  12)  using  an  array 
of  miniature  pressure  transducers.  These  measurements  were  carried 
out  in  a low  hub-to-tip  ratio  compressor  with  square-wave  and  some 
sine  wave  type  distortions.  A different  approach  was  taken  by 
Molkelke  (Ref.  13)  who  attempted  to  derive  the  unsteady  blade  lift 
from  measurements  of  absolute  flow  quantities.  He  investigated  the 
response  of  the  blading  to  square  wave  distortions  and  to  forward  and 
backward  ramp  type  distortions.  The  results  for  the  ramp  distortions 
were  found  to  differ  depending  on  which  way  the  total  pressure  defect 
was  penetrated.  The  unsteady  pressure  response  of  cascaded  blades  to 
gusts  was  examined  by  Satyanarayana  (Ref.  14)  in  a tunnel  where  the 
blades  remained  stationary  while  the  tunnel's  ceiling  and  floor  were 
moved.  The  resulting  gusts  were  primarily  non-convect ing  whereas  the 
inlet  distortion  problem  is  essentially  a convected  gust  problem. 

The  experimental  apparatus  also  restricted  the  investigation  to  very 
low  reduced  frequencies.  Kool , De  Ruyck,  and  Hirsch  (Ref.  15)  have 
made  measurements  of  the  wake  profiles  downstream  of  a compressor 
rotor  in  undistorted  flow  and  Colpin  and  Kool  (Ref.  16)  have  made 
similar  measurements  under  the  effects  of  a triangular  distortion 
pattern.  It  was  noted  that  the  wake  shapes  and  the  exit  flow  angles 
were  strongly  influenced  by  the  presence  of  the  distortion. 
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Some  of  the  work  which  has  been  done  on  compressor  instabilities 
will  now  be  considered.  As  was  previously  mentioned,  these  instabil- 
ities might  appear  as  surge  or  rotating  stall.  The  first  quantitative 
measure  of  when  a given  mode  will  occur  was  presented  by  Greitzcr 
(Ref.  2).  The  theoretical  derivation  which  leads  to  this  prediction, 
however,  requires  the  compressor  performance  curve  (Fig.  3)  as  an 
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FIG.  3 TYPICAL  COMPRESSOR  CHARACTERISTIC 

input.  Note  that  there  are  basically  two  branches  in  the  typical 
curve  shown.  These  correspond  to  unstalled  and  stalled  operation, 
respectively.  There  is  a range  of  inlet  conditions  where  either 


stalled  or  unstalled  operation  is  possible  with  a large  hysteresis 
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loop  apparent  between  stall  onset  (point  A to  point  B)  and  stall 
recovery  (point  C to  point  D) . Since  both  surge  and  rotating  stall 
are  strongly  dependent  on  the  shape  of  the  curve  in  the  vicinity  of 
this  hysteresis  loop,  it  is  important  to  know  the  overall  features  of 
the  performance  curve  for  both  stalled  and  unstalled  flow  including 
the  hysteresis  aspect.  A method  which  is  able  to  predict  the  basic 
features  of  these  hysteresis  loops  based  on  the  compressor  design 
parameters  is  described  in  Ref.  17.  At  present,  the  method  involves 
a great  deal  of  empiricism.  There  is  still  much  additional  work  to 
be  done . 

» 

The  rotating  stall  problem  has  also  been  the  subject  of  some 
study.  The  flow  in  a compressor  undergoing  rotating  stall  is  unsteady 
even  in  the  stationary  reference  frame.  Thus,  the  pneumatic  instrumen- 
tation commonly  used  to  measure  the  time-averaged  pressures  in  steady 
flows  is  not  generally  applicable.  As  pointed  out  by  Bennett  (Ref. 

18),  the  use  of  these  averaged  pressures  may  lead  to  significant 
errors  in  velocity  and  flow  angle  measurements  if  there  are  large 
variations  in  the  pressures  and  velocities  in  question.  The  diffi- 
culties involved  are  greatest  in  dealing  with  on-blade  measurements 
made  in  a non-uniform  flow  field.  Without  unsteady  instrumentation, 
one  can  only  obtain  average  pressures  (and  hence  normal  force  and 
pitching  moment)  at  an  average  inlet  angle.  This  type  of  data  does 


not  adequately  describe  the  aerodynamic  response  of  the  blade  to  the 
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distortion.  Pneumatic  instrumentation  can  be  used  to  measure 
performance  parameters  in  the  stationary  reference  frame  and  can  be 
quite  useful  when  used  in  conjunction  with  high  response  instrumenta- 
tion (c.f.  the  experimental  work  reported  in  Ref.  2).  Sexton, 

O'Brien  and  Moses  (Ref.  19)  have  developed  a telemetry  system  and 
have  used  it  in  conjunction  with  miniature  pressure  transducers  to 
make  on-rotor  studies  of  rotating  stall.  This  system  is  necessarily 
channel  limited  and  has  a restricted  bandwidth. 

It  is  obvious  that  further  extension  of  the  analytical  efforts 
previously  described  requires  the  availability  of  good  experimental 
data  both  as  an  input  to  the  theories  and  as  a check  on  their  accuracy. 
Therefore,  North  Carolina  State  University  (NCSU)  and  United 
Technologies  Research  Center  (UTRC)  have  embarked  on  a joint  project 
under  the  sponsorship  of  the  Air  Force  Office  of  Scientific  Research 
( AFOSR)  to  study  analytically  and  experimentally  the  response  of  a 
compressor  stage  to  an  inlet  distortion.  The  work  reported  herein  is 
the  first  of  several  planned  experimental  investigations  of  this 
phenomenon  using  UTRC's  Large  Scale  Rotating  Rig  (LSRR). 


EXPERIMENTAL  FACILITY 


Large  Scale  Rotating  Rig 

This  experiment  was  conducted  in  the  Large  Scale  Rotating  Rig 
(LSRR)  located  at  United  Technologies  Research  Center  (UTRC)  in  East 
Hartford,  Connecticut.  The  LSRR  was  originally  designed  to  provide  a 
facility  which  would  be  of  sufficient  size  to  permit  resolution  of 
endwall  flows,  possess  a high  degree  of  flexibility  in  regard  to  the 
configurations  which  can  be  tested,  and  enable  measurements  to  be 
made  directly  in  the  rotating  reference  frame. 

Overall  arrangement  of  the  facility  is  shown  in  Fig.  4.  Flow 
enters  the  rig  through  a twelve  foot  diameter  inlet.  A six  inch 
thick  section  of  honeycomb  material  is  mounted  at  the  inlet  face  to 
remove  any  crossflow  effects.  The  inlet  contracts  smoothly  down  to  a 
cross  section  five  feet  in  diameter  with  the  flow  then  passing 
through  a series  of  fine  mesh  screens.  Immediately  downstream  of  the 
screens  is  a telescoping  section  which  slides  axially  to  permit 
access  to  the  test  section.  The  test  section  is  five  feet  in  diameter 
with  large  segments  of  the  outer  casing  constructed  of  Plexiglas  to 
facilitate  flow  visualization,  holography,  or  laser-velocimeter 
studies.  Axial  length  of  the  test  section  (excluding  the  hub  spinner 
length)  is  three  feet.  The  rotor  is  mounted  to  a shaft  cantilevered 
from  two  downstream  bearings  thus  providing  a clean  flow  path  to  the 
most  upstream  blade  row.  A hydraulic  motor  capable  of  turning  the 
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FIG.  4 UTRC  LARGE  SCALE  ROTATING  RIG 
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shaft  at  800  rpm  drives  the  rotor.  Downstream  of  the  test  section, 
flow  passes  through  an  annular  diffuser  into  a centrifugal  fan  and  is 
subsequently  exhausted  from  the  rig.  A vortex  valve  mounted  at  the 
fan  inlet  controls  the  flow  rate. 

The  configuration  tested  in  this  experiment  was  an  isolated 
rotor  with  a hub-to-tip  ratio  of  0.8.  The  rig  had  previously  been 
used  only  with  models  which  involved  a stationary  set  of  inlet  guide 
vanes  upstream  of  the  rotor  and  these  vanes  had  been  used  to  support 
the  hub  spinner  which  was  also  stationary.  Furthermore,  no  provision 
existed  for  mounting  a distortion  generating  screen.  The  modifica- 
tions to  the  rig  necessary  to  support  the  hub  spinner  without  the 
inlet  guide  vanes  and  to  provide  a place  for  mounting  the  distortion 
screen  are  shown  in  Fig.  5.  A large  bearing  was  added  on  the  front 
of  the  rotor  shaft  to  support  the  rear  of  the  hub  spinner.  A spacer 
ring  was  installed  between  two  segments  of  the  hub  spinner  so  the 
forward' junct ion  coincided  with  the  joint  in  the  outer  casing  between 
the  test  section  and  the  telescoping  section  upstream.  The  distortion 
screen  was  installed  at  this  location  to  support  the  front  of  the 
spinner  and  restrain  it  from  rotating.  When  the  rig  was  run  without 
distortion,  a screen  frame  (without  any  wire  screening  installed)  was 
used  in  place  of  the  distortion  screen  to  perform  this  support 


function . 
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ROTATING  MEMBERS 

PARTS  ADDED  FOR  INLET 
DISTORTION  EXPERIMENT 


FIG.  5 MODIFICATIONS  TO  LSRR  TO  PERMIT  INSTALLATION  OF 
DISTORTION  SCREENS 
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The  rotor  had  28  blades  with  a span  and  chord  of  0.5  feet.  Each 
blade  consisted  of  an  eight  percent  NACA  65-series  thickness  distribu- 
tion superimposed  on  a 50-degree  circular-arc  mean  camber  line.  The 
blades  were  installed  with  a midspan  stagger  angle,  B*,  of  30  degrees. 
The  gap-to-chord  ratio  at  midspan  was  1.01.  The  blade  cross  section 
at  midspan  and  the  relative  frame  velocity  and  flow  angle  nomenclature 
are  shown  in  Fig.  6. 

I 

The  LSRR  was  operated  at  a constant  corrected  wheel  speed  of  600 
rpm  throughout  the  experiment.  The  mean  :nle:  angle  was  varied  by 
changing  the  flow  rate  with  the  downstream  vortex  valve.  The  circum- 
ferentially averaged  axial  velocity  in  the  test  section  was  varied 
between  75  ft/sec  and  170  ft/sec,  resulting  in  a blade  Reynolds 
number  of  approximately  0.5  x 10^.  This  value  is  comparable  to 
that  found  in  an  actual  turbomachine.  Therefore,  although  the  flow 
is  admittedly  incompressible,  both  the  blade  geometry  and  the  Reynolds 
number  are  representative  of  the  blading  in  the  high  compressors  of 
current  engines. 

Distortion  Screens 

The  inlet  distortion  imposed  on  the  rotor  in  this  experiment  was 
a nominal  20-percent,  once-per-revolution,  sinusoidal  variation  in 
the  inlet  axial  velocity.  This  distortion  was  generated  by  a screen 
consisting  of  segments  of  wire  screening  of  various  meshes  attached 
to  a support  frame.  i- " ssign  of  the  distortion  screen  follows  the 
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work  of  Bruce  (Ref.  20).  Using  the  equations  presented  in  Ref.  20,  a 
resistance  coefficient  distribution  for  the  screen  was  calculated 
which  would  give  the  desired  distortion.  Various  meshes  of  wire 
screening  were  chosen  so  that  twelve  different  values  of  resistance 
coefficient  could  be  obtained  using  one  or  more  of  the  meshes.  The 
desired  distribution  and  that  actually  obtained  are  shown  in  Fig.  7. 

Instrumentation  System 

One  of  the  most  important  aspects  of  any  experimental  program  is 
the  instrumentation  system.  It  must  provide  accurate  measurements 
with  a minimum  of  attention.  Part  of  the  instrumentation  used  in 
this  experiment  was  mounted  inside  the  rotor  hub  of  the  LSRR  and 
rotated  with  the  rotor,  thus  imposing  additional  constraints  on  size 
and  durability.  Furthermore,  the  limited  number  of  slip  rings 
available  to  interface  the  rotating  equipment  with  the  stationary 
equipment  was  a key  factor  in  many  of  the  decisions  made  during  the 
design  phase.  The  system  described  below  was  the  result  of  consider- 
able experimentation.  While  no  single  component  was  totally  new,  the 
manner  in  which  the  various  pieces  of  equipment  were  combined  was 
unique  and  effective.  Some  of  the  techniques  originally  developed 
for  the  inlet  distortion  experiment  have  found  application  in  other 


experimental  programs  at  UTRC. 
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Pneumatic  Instrumentation 


Normal  instrumentation  for  the  LSRR  consists  of  stationary  and 
rotating  reference  frame  pneumatic  probes  and  taps  connected  to 
scani-valves  to  measure  the  time  averaged  pressures  at  various  points 

in  tlie*  flow  field.  The  stationary  trame  pneumatic  lnStlUlUHiilHLiuii 

was  ^mplcyed  to  measure  static  and  total  pressurM4-*n  an  di^ht  point 
<^rcumfei  ential  array  upstream  of  the  rotor  and^the  static  prrmsure 
pi  a ^imi  lar  array  downstream  of  the  rotor,  7*low  angles  upstrefflfrj and 

m A.  . . \ 

ftwnstrei  m of  the  rotor  were  measured  at  afTsingie  circumferential^ 

« / A 

location  using  pneumatic  wedge  probes,  /(The  location  of  the  stationary 

n 1/  " 

v>  \ 

Jfiram26in  trumentation  is  illustrated  jtn  Fig.  8.  In  the  rotating 
a.  A DESIRED  DISTRIBUTION 

refe^nc  —fAoe,  one  blade  was  ins^menAd^f^^p^i^^^l®^  in  a 
Gaussian  array  vtef.  21)  at  i ts/mi dspan . The  details  of  the  six- 


element  array  emplo> 


iis  experiment  and  the  measuring  station 


nome  nyla^  ure  is|  given  |in  Fig^  9 and  | Table  J. , By  prescribing  tl|e  j 

station^  at  wh46h  the*t>ressdWs  onlB&ie  suJQffece  of4fthe  tSBBde  a 880  to  tJ0O 

measured,  large  gains  may  be  realized  in  the  accuracy  of  the  inte- 

CIRCUMFERENTIAL  ANGLE  ,0.  DEG 

grated  normal  force  and  pitching  moment  coefficients.  An  accuracy  on 
the  order  of  three  percent  should  be  realized  with  the  six-element 

arFW.  ^efiBSIsnh*NeB>edeFPteto9t  D?#tRIBUTION  OF  DISTORTION  SCREEN 

Unsteady  Instrumentation 

To  measure  the  fluid  dynamic  response  of  the  rotor  to  the  inlet 
distortion,  it  was  necessary  to  provide  high  response  instrumentation 


TABLE  I 


GAUSSIAN  POSITIONS  AND  WEIGHTING  FACTORS 


Stat ion 

1,7 

0.0337652 

0.0856622 

2,8 

0.1693953 

0.1803808 

3,9 

0.3806904 

0.2339570 

4,10 

0.6193096 

0.2339570 

5,11 

0.8306047 

0.1803808 

6,12 

0.9662348 

0.0856622 

Note : 

Station  numbers  1-6  are 

on  the  suction  surface 

and  numbers  7-12  are  on 

the  pressure  surface. 
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in  the  rotating  reference  frame.  Since  no  such  instrumentation 
existed,  it  was  necessary  to  design  a package  which  could  make  the 
desired  measurements.  Parameters  of  interest  were: 

1.  Rotor  blade  surface  pressures 

2.  Rotor  blade  transient  surface  flow  phenomena 

3.  Pressures  and  velocities  in  the  near  wake  of  the  rotor. 

The  rotating  frame  instrumentation  in  the  LSRR  was  interfaced  to 
the  stationary  frame  through  a slip  ring  package.  A system  was 
evolved  which  maximized  the  number  of  channels  of  information  which 
could  be  obtained  with  the  limited  number  of  slip  rings  available  and 
simultaneously  minimized  the  noise  contamination  problems  commonly 
experienced  with  slip  rings.  The  basic  idea  was  to  place  all  the 
signal  conditioning  and  amplifying  equipment  necessary  to  produce  a 
high  level,  single  ended  signal  inside  the  rotor  hub.  After  a 
limited  number  of  slip  rings  were  dedicated  to  supplying  power  and  a 
common  output  reference,  only  one  slip  ring  was  needed  to  transmit 
each  channel  of  data.  Furthermore,  by  keeping  the  current  through 
the  slip  rings  which  transmit  the  data  to  a minimum,  it  was  possible 
to  virtually  eliminate  the  electrical  noise  normally  generated  by  the 
varying  contact  resistance  which  is  inherent  in  the  slip  rings. 

Power  Supplies 

An  unregulated  direct  current  power  supply  lcoated  in  the 
stationary  frame  produced  positive  and  negative  outputs  of  twenty 
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vo 1 1 8 with  respect  to  a common  ground  with  a current  capability  of 
three  amperes.  The  power  was  brought  onto  the  rotating  frame  through 
the  slip  rings  where  it  was  regulated  to  produce  the  necessary 
voltages  for  driving  the  various  pieces  of  equipment.  As  shown  in 
Fig.  10,  three  regulators  which  produce  plus  and  minus  fifteen  volts 
were  used  to  power  the  constant  temperature  anemometers  and  pressure 
transducer  amplifiers.  A separate  ten  volt  regulator  provided 
excitation  to  the  pressure  transducers  themselves.  Schematics  for 
the  regulators  are  presented  in  Appendix  I. 

Surface  Pressure  Measuring  Systems 

Unsteady  pressures  on  the  surface  of  one  of  the  rotor  blades 
were  measured  using  miniature  pressure  transducers.  The  units  used 
in  this  experiment  were  Kulite  LQ-64-125  + 5D  transducers  which  have 
an  exceptionally  high  sensitivity  at  the  expense  of  a slightly  higher 
than  normal  drift  rate.  The  installation  of  a typical  transducer  is 
shown  in  Fig.  11.  The  transducers  which  measure  the  pressures  on  the 
suction  surface  of  the  blade  were  inserted  from  the  pressure  surface, 
and  vice  versa.  The  transducers  were  installed  with  Duco  cement  and 
pressure  sealed  with  model  airplane  dope.  The  depressions  in  which 
the  transducers  rest  were  filled  with  a mixture  of  model  airplane 
dope  and  talcum  powder  almost  to  the  surface  of  the  airfoil.  When 
this  compound  had  cured  sufficiently,  a layer  of  epoxy  putty  was 


applied  and  sanded  to  the  airfoil  contour.  This  process  resulted  in 
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FIG.  11  INSTALLATION  OF  PRESSURE  TRANSDUCER  IN  BLADE 

an  extremely  smooth  surface  which  offered  minimal  disturbance  to  the 
flow.  Also,  transducers  were  readily  recovered  by  chipping  away  the 
epoxy  layer  and  dissolving  the  remainder  of  the  potting  material  with 
acetone . 

As  was  previously  mentioned,  the  regulator  supplying  excitation 
to  the  transducer  was  mounted  inside  the  rotor  hub.  Each  transducer 
was  connected  to  a thermal  compensation  network  (Fig.  12)  which  was 
designed  to  minimize  thermal  induced  drift.  The  transducer  output 
was  connected  to  an  amplifier  with  a fixed  gain  of  100.  The  circuit 
for  this  amplifier  is  included  in  Appendix  II  along  with  a brief 
circuit  description.  The  amplified  signal  was  brought  off  the 
rotating  frame  through  the  slip  ring  assembly  to  the  stationary  frame 
instrumentation.  The  time  averaged  portion  of  the  signal  (which  was 
not  needed  because  the  time  averaged  pressures  were  measured  with 
pneumatic  instrumentation)  were  removed  by  passing  the  signal  through 
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FIG.  12  TYPICAL  CIRCUITRY  FOR  PRESSURE  TRANSDUCER  CHANNELS 
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a high-pass  filter  with  an  appropriate  cutoff  frequency.  A calibra- 
tion relay  allowed  either  this  filtered  signal  or  a calibration 
voltage  to  be  applied  to  the  input  of  a variable  gain  instrumentation 
amplifier  so  that  the  system  could  be  recalibrated  every  time  the 
instrumentation  amplifier  gain  was  changed.  (The  calibration  proce- 
dures will  be  discussed  in  more  detail  in  a subsequent  section.)  The 
output  of  each  of  these  instrumentation  amplifiers  was  connected  to 
one  of  the  tape  recorder  inputs. 

Skin  Friction  Gages 

Transient  surface  flow  phenomena  were  monitored  by  Micro- 
Measurements  EGT-50  heated-film  skin  friction  gages  bonded  to  the 
surface  of  the  blade  adjacent  to  each  of  the  pressure  transducer 
locations.  Wiring  for  the  skin  friction  gages  was  kept  below  the 
surface  of  the  blade  to  minimize  flow  disturbances  due  to  the  instru- 
mentation. The  gages  operated  in  the  constant  temperature  (constant 
resistance)  mode  by  anemometer  modules  located  inside  the  rotor  hub 
(Fig.  13).  These  anemometers  were  designed  especially  for  this 
experiment  and  were  made  as  simple  and  durable  as  possible  so  that 
they  could  withstand  the  acceleration  loads  encountered  in  the  rotor 
hub.  A schematic  diagram  and  brief  description  of  the  circuit  are 
included  in  Appendix  III.  The  output  of  an  anemometer  module  was 
relatively  high  in  amplitude  and  needed  no  further  amplification. 


After  the  signal  had  passed  through  the  slip  rings,  the  time  averaged 
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component  was  removed  by  a low  pass  filter.  No  attempt  was  made  to 
recover  this  time  averaged  part  because  the  skin  friction  gages  were 
used  only  for  qualitative  information.  The  unsteady  portion  of  the 
signal  was  then  routed  through  the  calibration  relay  to  the  recording 
equipment . 

Rotating  Frame  Wake  Probe 

A wake  probe  consisting  of  a triaxial  hot-film  probe  to  measure 
velocity  and  flow  angles  and  a pressure  transducer  to  measure  total 
pressure  was  mounted  immediately  behind  the  rotor  in  the  rotating 
reference  frame  (Fig.  14).  The  hot-film  probe  was  a modified  unit 

TOTAL  PRESSURE  PROBE 

TRIAXIAL  HOT -FILM  PROBE 


FIG.  14  WAKE  PROBE  USED  IN  LSRR 

manufactured  by  Thermo  Systems,  Inc.  and  designated  TSI  1296  M.  Its 
sensors  were  operated  in  the  constant  temperature  mode  by  anemometers 
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of  the  same  type  as  those  used  with  the  skin  friction  gages  (Appendix 
III).  These  anemometers  were  packaged  differently  because  they  were 
mounted  in  a different  area  of  the  hub  than  the  others.  As  in  the 
case  of  the  skin  friction  gages,  the  anemometer  output  needed  no 
further  amplification  and  was  routed  through  the  slip  rings  to  the 
stationary  frame  instrumentation.  The  time  averaged  portions  of  the 
signals  were  read  from  a digital  voltmeter  and  manually  recorded. 

Thus,  even  though  the  time-averaged  portions  were  not  recorded  on 
tape,  they  were  available  so  that  the  complete  signals  could  be 
subsequently  reconstructed.  The  remainder  of  the  circuitry  was 
identical  to  that  used  for  the  skin  friction  gages. 

The  total  pressure  in  the  wake  was  measured  by  a Kulite  XCQL-093- 
25A  pressure  transducer  mounted  in  a sharp  edged  tube  (Fig.  15).  The 


FIG.  15  KULITE  TOTAL  PRESSURE  PROBE 
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circuitry  used  with  this  transducer  was  identical  to  that  employed  with 
the  blade  pressure  transducers.  In  this  particular  instance,  it  was 
desired  to  obtain  the  steady  as  well  as  the  unsteady  portion  of  the 
total  pressure  from  the  transducer;  consequently,  the  high-pass  filter 
should  have  been  switched  out  of  the  circuit.  However,  it  was  acci- 
dentally left  in  the  circuit  and  only  the  unsteady  part  of  the  signal 
was  recorded . 

Data  System 

Unsteady  data  were  recorded  on  the  UTRC  Analog  Recording  System 
(ARES),  consisting  of  two  fourteen-channel  FM  tape  recording  units  and 
their  associated  amplifying  and  control  equipment.  (The  calibration 
relays,  instrumentation  amplifiers,  and  tape  units  indicated  in  Figs. 

12  and  13  wfere  actually  part  of  ARES.)  A modified  IKIG  B time  code  was 
recorded  on  one  channel  of  each  tape  unit  to  permit  time  correlation  of 
the  two  tapes  and  to  facilitate  record  search  and  playback  during 
subsequent  data  processing.  The  use  of  this  time  code  allows  the  data 
from  the  two  tapes  to  be  synchronized  to  within  10-^  seconds.  The 
analog  data  contained  on  the  ARES  tapes  were  digitized  and  recorded  on 
digital  magnetic  tape  prior  to  data  reduction. 

The  pneumatic  data  and  test  conditions  were  measured  by  the  Test 
Stand  42  static  data  system.  This  system  writes  the  data  directly  onto 
digital  magnetic  tape.  The  data  are  then  processed  using  built-in  cal- 
ibration functions  to  provide  the  various  parameters  in  physical  units. 
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SCOPE  OF  EXPERIMENT  1 

1 

| 

To  place  the  present  experiment  in  its  proper  perspective,  it  is 
useful  to  examine  the  possible  modes  in  which  a compressor  may  ' 

operate  using  Fig.  16.  There  are  four  main  regimes  in  which  operation 
is  possible: 


FIG.  16  FLOW  REGIMES  FOR  COMPRESSOR  OPERATION 

I)  undistorted  flow  with  no  rotating  stall 

II)  distorted  flow  with  no  rotating  stall 

III)  undistorted  flow  with  rotating  stall  occuring 

IV)  distorted  flow  with  rotating  stall  occuring 
There  are  also  two  transition  regions  (l-II  and  II-IV)  in  which  the 
compressor  may  move  in  and  out  of  rotating  stall.  Region  I is  the 
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design  operating  condition  of  the  compressor.  The  previously  discussed 
problems  may  be  seen  to  fall  in  regions  II,  III,  IV  and  the  transition 
region  II-IV.  Of  these  regions,  IV  is  the  one  in  which  measurements 
are  least  meaningful,  extremely  difficult,  and  subject  to  a variety 
of  errors.  Even  if  measurements  are  made  accurately,  it  is  doubtful 
that  the  distortion-induced  stall  mode  measured  in  a test  rig  will 
accurately  represent  that  in  an  actual  engine.  This  is  because  the 
distortion  which  induced  the  engine  stall  will  often  be  subsequently 
removed  so  that  the  compressor  is  stalled  even  though  the  inlet  flow 
is  now  uniform.  The  most  suitable  topics  for  investigation  are  thus 
compressor  operation  in  rotating  stall  with  no  imposed  distortion 
(region  III)  and  compressor  operation  in  distorted  flow  but  without 
rotating  stall  (region  II).  This  latter  topic  may  be  extended  to 
include  the  transition  region  II-IV.  The  capabilities  of  the  LSRR 
(and  the  analytical  model  being  developed  concurrently  at  NCSU 
(unpublished))  are  best  suited  for  investigation  of  both  of  these 
regions  (II  and  II-IV).  During  this  first  experimental  phase, 
measurements  were  made  in  all  of  the  regions  in  order  that  the 
operating  limitations  of  the  LSRR  and  the  capability  of  the  available 
instrumentation  to  make  meaningful  measurements  could  be  evaluated. 
However,  following  the  reasoning  outlined  above,  measurements  were 
concentrated  in  region  1 to  provide  a baseline  case  and  in  region  II 


to  obtain  unsteady  data. 
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It  is  well  known  that  an  airfoil  in  a subsonic  stream  influences 
the  upstream  flow  as  well  as  the  downstream  flow.  Extending  this 
concept  to  a compressor,  it  is  clear  that  the  interaction  between  the 
rotor  and  the  incoming  distorted  flow  field  establishes  a velocity 
distribution  at  the  compressor  face  that  determines  the  inlet  flow  to 
which  individual  compressor  blades  respond.  From  a gap-averaged 
point  of  view,  the  unsteady  response  of  a rotor  blade  row  is  dependent 
on  the  imposed  inlet  angle  variations.  Thus,  a rotating  blade  in  a 
multi-stage  machine  may  be  thought  of  as  being  affected  by  its 
neighboring  blade  rows  primarily  through  the  imposed  or  induced 
upstream  flow  angle  distribution.  The  measurement  and  understanding 
of  the  variations  in  the  performance  of  a single  blade  row  with  inlet 
angle  is  therefore  a necessary  prelude  to  the  understanding  of  the 
complex  interactions  in  a multistage  machine.  The  performance  of  a 
blade  row  may  be  quantified  in  terms  of  the  variations  in  flow  exit 
angle  and  total  pressure  loss  with  changing  inlet  angle.  This  inlet 
angle  distribution  may  be  simply  the  result  of  an  imposed  distortion 
or  it  may  arise  from  the  interactions  of  the  rotor  with  other  compo- 
nents such  as  downstream  stators,  diffusers,  struts,  etc.  Of  course, 
the  blade  row  performance  may  be  more  easily  and  more  completely 
understood  if  the  inlet  angle  variation  is  simple.  In  the  past, 
various  investigators  have  studied  sinusoidal  distortions,  aquare-wave 
distortions  and  ramp  or  triangle-wave  distortions.  The  present 
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experimental  program  was  aimed  at  measuring  the  response  of  an 
isolated  rotor  to  a sinusoidal  variation  in  inlet  axial  velocity. 

This  choice  was  made  for  a number  of  reasons.  First,  the  sinusoidal 
distortion  is  relatively  easy  to  produce  (Ref.  20).  Second,  the 
square-  and  triangular-wave  distortions  may  be  represented  in 
terms  of  the  sine-wave  by  means  of  a Fourier  Series.  The  study  of 
sinusoidal  distortion  patterns  of  various  reduced  frequencies  will 
thus  yield  information  pertinent  to  other  types  of  distortions. 
Finally,  the  square-  and  triangular-wave  distortions  may  be  thought 
of  as  having  higher  order  harmonics  which  may  make  the  understanding 
of  the  rotor  response  more  difficult.  Conversely,  the  sine-wave 
distortion  contains  only  one  harmonic.  The  reduced  frequency  may  be 
readily  controlled  by  varying  the  number  of  cycles  of  the  distortion 
pattern  the  screen  produces.  In  the  present  experiment,  a one-cyc le- 
per-revolution screen  (and  the  attendant  low  reduced  frequency)  was 
chosen  to  produce  some  unsteady  effects  without  causing  a significant 
departure  from  the  steady  state  behavior. 
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EXPERIMENTAL  PROCEDURE 
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Data  Acquisition 

During  this  experiment,  the  LSRR  was  operated  at  a constant 
corrected  wheel  speed  of  600  rpm.  For  a given  upstream  total  tem- 
perature, Tq  , the  actual  rotational  speed  required  to  achieve 
the  corrected  wheel  speed  of  600  rpm  is  given  by  Eq.  (1). 

T 

Nact  = 600  (1) 

518. 7°R 

The  rotor  rpm  was  measured  by  a digital  tachometer  and  was  found 
to  be  stable  once  it  had  been  set.  The  axial  velocity  in  the  case  of 
undistorted  flow  was  varied  so  that  the  blade  incidence  angle  (refer- 
enced to  the  mean  camber  line)  ranged  from  a high  negative  value 
(a^i  * -18*)  up  to  the  point  of  rotating  stall  inception  ( amc  ^ = 1.5*). 
In  distorted  flow,  the  mean  incidence  angle  was  varied  over  approxi- 
mately the  same  range.  Due  to  the  20  percent  distoriton  in  inlet 
axial  velocity,  the  instantaneous  inlet  angle  varied  on  the  order  of 
♦ 6*  about  the  mean  angle.  In  all  cases  tested,  the  stall  mode  was 
observed  to  be  a multiple-cell,  partial-span  rotating  stall  near 
the  rotor  hub.  To  set  the  desired  axial  velocity,  a U-tube  manometer 
was  connected  between  a wall  static  tap  and  a total  pressure  probe 
located  in  the  center  of  the  annulus.  The  axial  position  of  both  the 
static  tap  and  the  total  pressure  probe  was  midway  between  the  rotor 


and  the  distortion  screen.  The  measuring  station  was  located  at  the 
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90°  point  in  the  distortion  pattern  where  the  local  axial  velocity 
would  be  equal  to  the  average  axial  velocity  if  there  were  no  preswirl. 

A more  accurate  measurement  of  the  circumferential  distribution  of 
axial  velocity  was  made  using  a circumferential  array  of  measuring 
stations.  This  axial  velocity  distribution  was  used  in  reducing  the 
data,  but  the  visual  reading  obtained  from  the  single  U-tube  manometer 
proved  adequate  for  setting  the  operating  point. 

The  number  of  channels  which  could  be  simultaneously  transmitted 
from  the  rotating  frame  was  limited  to  17  by  the  slip  ring  assembly. 

It  was  thus  necessary  to  run  two  instrumentation  configurations  to 
obtain  the  desired  data.  The  first  configuration  obtained  data  from 
the  complete  set  of  pressure  transducers  plus  selected  skin-friction 
gages.  The  pressures  were  integrated  to  obtain  blade  normal  force 
and  pitching  moment  coefficients.  The  second  instrumentation  confi- 
guration included  all  the  suction  surface  pressure  transducers  and 
skin  friction  gages  as  well  as  the  wake  probe.  For  a given  flow 
condition,  the  wake  probe  was  traversed  across  about  1.5  blade  gaps 
giving  particular  emphasis  to  the  blade  wake  regions.  Due  to  the 
time  involved  in  making  this  traverse  and  the  mass  of  data  which 
would  be  obtained,  only  a few  flow  conditions  were  run  in  this 

configuration.  j 

The  instrumentation  available  in  each  of  the  two  configurations 

j 

is  presented  in  Appendix  IV  which  also  includes  a chronological 

I 

I 


I 

I 
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listing  of  the  test  points  run. 


Notes  are  included  to  point  out 
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problems  encountered  during  the  tests  which  may  influence  the  accuracy 
of  the  results. 

Once  a desired  operating  condition  had  been  reached  and  the  flow 
had  stabilized,  the  gain  of  each  channel  of  the  unsteady  data  system 
was  set  as  high  as  possible  without  overloading  the  tape  recording 
equipment.  An  oscilloscope  was  used  to  monitor  the  tape  inputs  to 
insure  that  the  signal  was  within  the  required  limits.  Calibration 
records  were  then  acquired  for  both  the  static  and  unsteady  data 
systems  and  the  data  were  recorded  immediately  thereafter.  If  a wake 
traverse  were  being  made  with  the  rotating  frame  wake  probe,  static 
data  would  be  taken  only  at  the  first  point  for  each  flow  condition. 

The  data  from  the  static  system  was  subsequently  converted  to 
pressures  and  temperatures  by  a computerized  process  making  use  of 
the  calibration  records.  The  unsteady  data  were  digitized  by  the 
process  of  Ref.  23.  After  completion  of  the  digitizing  process,  the 
calibration  records  were  used  to  convert  the  data  to  appropriate 
physical  units. 

Overall  Performance  Parameters 

The  wall  static  and  mid-span  total  pressures  were  obtained  at 
eight  circumferential  stations  45'  apart  approximately  1.5  chord 
lengths  upstream  of  the  rotor.  Using  the  Bernoulli  equation,  the 


velocity  at  each  of  these  was  calculated.  The  velocity  was  assumed 
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to  be  purely  axial  at  these  locations.  Limited  flow  angle  measure- 
ments in  undistorted  flow  did  verify  this  assumption.  These  veloci- 
ties were  then  fitted  to  the  form  given  below  in  Eq . (2)  by  a least 
squares  technique. 


Cx  (6)  = Cx  + A cos  (6  + 4>) 


(2) 


The  angle  0 is  the  circumferential  position  in  the  distortion  pattern 
and  increases  in  the  direction  of  rotation  of  the  rotor.  The  quantity 
Cx  is  the  average  value  of  the  axial  velocity  and  the  constants  A and 
C are  the  amplitude  and  phase  angle  of  the  velocity  variation,  respec- 
tively. The  static  pressure  upstream  of  the  rotor  was  fitted  to  the 
form  given  by  Eq . (3). 


P (0)  = P + B COS  (6  + <J»)  , 

Si  Si 


(3) 


where  P is  the  average  value  of  the  static  pressure  and  B and  i|<  are 
S1 

the  amplitude  and  phase  angle  of  the  pressure  variation.  The  Bernoulli 
equation  was  applied  once  again  to  obtain  a smoothed  total  pressure 
profile  around  the  distortion  which  is  consistent  with  the  velocity 
and  static  pressure  profiles  obtained  from  the  curve  fits.  A small 
pressure  and  velocity  distortion  was  found  to  exist  even  with  the 
distortion  screen  removed.  The  velocity  distortion  without  the 
screen  varied  between  two  and  four  percent.  With  the  screen  in 
place,  this  pre-existing  distortion  was  found  to  interact  with  that 
of  the  screen  to  produce  a net  distortion  of  lower  amplitude  than  had 
been  anticipated.  For  all  of  the  results  presented,  the  measured 
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upstream  conditions  as  smoothed  by  the  previously  discussed  techniques 
are  used.  The  inlet  flow  angle,  Bj  , is  based  on  the  axial  velocity 
at  the  upstream  measuring  station  and  does  not  include  any  swirl  that 
may  have  occurred  between  that  axial  position  and  the  rotor  inlet 
face.  The  inlet  angle  is  given  by 


6 i = tan 


(4) 


For  distorted  inlet  flow,  the  inlet  angle  will  vary  with  position  in 
the  distortion  pattern;  however,  an  average  inlet  angle,  may  be 

calculated  by  inserting  the  average  axial  velocity,  Cx , into  Eq.  (4). 

The  exit  angle,  62,  is  an  important  measure  of  the  performance 
of  the  rotor.  For  undisturbed  flow,  this  angle  was  computed  from  the 
measured  stationary  frame  flow  angle  downstream  of  the  rotor,  . 

The  computation  was  iterative  and  began  by  assuming  that  the  axial 
velocity  downstream  of  the  rotor  was  identical  to  that  upstream.  The 
exit  angle  is  given  by 


&2  = tan 


-1 


x2 


U - CX2  tan©2 


(5) 


where  C is  the  downstream  axial  velocity  at  midspan.  Note  that 
x2 

the  spanwise  averaged  axial  velocities  upstream  and  downstream  of  the 
rotor  must  be  identical  to  satisfy  conservation  of  mass.  However, 
due  to  the  thickening  of  the  wall  boundary  layers  at  the  root  and  tip 
of  the  blading,  the  midspan  axial  velocity  downstream  of  the  rotor 


will,  in  general,  be  greater  than  that  upstream  of  the  rotor.  The 
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importance  of  this  phenomenon  was  not  fully  realized  at  the  outset  of 
the  experiment  and  no  provision  was  made  for  measuring  the  axial 
velocity  downstream  of  the  rotor.  Instead,  it  was  assumed  to  be  the 
same  as  that  upstream  of  the  rotor.  When  the  magnitude  of  the  error 
thus  introduced  was  realized,  it  became  necessary  to  devise  a method 
of  estimating  the  downstream  axial  velocity.  Horlock  (Ref.  24) 
presents  a correlation  between  the  axial  velocity  ratio  across  a 
cascade  and  the  flow  turning  due  to  Rhoden  (Ref.  25).  Using  this 
correlation,  one  may  iterate  to  obtain  the  correct  exit  angle.  As  a 
check  on  the  validity  of  this  correction,  the  exit  angle  was  also 
calculated  using  a simple  momentum  theory  (see  Appendix  V). 

Another  important  performance  parameter  is  the  rotor  pressure 
rise  coefficient  which  is  defined  as 

p s 2 - pS! 

pr  = j (6) 

1 / 2 pU 

where  P is  the  downstream  static  pressure,  p is  the  fluid 
s2 

density,  and  U is  the  tangential  wheel  speed.  To  be  consistent,  this 
coefficient  must  be  defined  at  the  blade  midspan.  Therefore,  the 
downstream  static  pressure  must  be  corrected  (from  a wall  value) 
using  the  radial  equilibrium  equation. 


if,  - p cfl2 

dr  r 


(7) 
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Integrating,  we  obtain 


tip 


-p  -/’ 
1 H 


'“V2  dr 
pc0  ~ , 


(8) 


where  i®  the  pressure  at  the  blade  tip  (outer  casing)  and  Pj 

is  the  pressure  at  the  radial  station  in  question.  The  radii,  r^ 
and  rtip,  are  those  of  the  station  in  question  and  the  tip,  respec- 
tively. This  equation  may  be  greatly  simplified  by  assuming  p and 


Cq  to  be  constants 


ptiP  - pi  “ pce  ln(“^£) 


(9) 


and  the  result  may  be  used  to  calculate  the  static  pressure  at  the 
blade  midspan  from  measured  pressures  on  the  outer  casing. 

A final  important  parameter  used  to  measure  the  performance  of  a 
compressor  stage  is  the  loss  coefficient.  For  undistorted  flow,  the 
steady-state  loss  coefficient  is  given  by 


Xss 


(10) 


where  P’  and  P’  are  the  total  pressures  upstream  and  down- 
°1  °2 

stream  of  the  rotor,  respectively,  in  the  rotating  reference  frame 
and  Wj  is  the  relative  velocity  upstream  of  the  blading  (see 
Fig.  6). 


Unsteady  Blade  Pressures 

Reduction  of  the  unsteady  data  requires  that  the  pressures 
measured  on  the  surface  of  the  blade  be  converted  to  an  appropriate 
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nondimen8ional  form.  The  most  relevant  parameter  to  use  is  the 
pressure  coefficient  defined  by 


C = 


P - P. 


*5PW  j 


(11) 


where  P is  the  static  pressure  upstream  of  the  blade  leading 
81 

edge.  With  distorted  inflow,  this  parameter  varies  as  the  blade 
traverses  the  distortion.  The  relative  velocity,  Wj , is  used  to 
nondimens ionalize  the  pressure  coefficients.  It  should  be  noted  that 
is  not  constant  but  varies  around  the  distortion.  The  computed 
coefficients  are  phase-lock  averaged  over  five  cycles  (Ref.  26)  to 
minimize  the  effects  of  any  random  fluctuations  in  the  pressure  data. 
A Fourier  series  representation  of  the  average  pressure  coefficient 
variation  at  each  chordwise  station  was  computed  using  the  formula 


Cp  (6) 


N 

a0  + I [an  cos  n 9 + bn  sin  n 0] 
n=l 


(12) 


and  the  coefficients  of  this  series  were  then  filtered  with  a cosine- 
square,  low-pass  filter  (Ref.  27).  The  cutoff  frequency  of  the 
filter  was  chosen  to  retain  the  response  of  the  blade  pressures  to 
the  distortion  but  to  remove  the  response  of  these  pressures  to 
screen  irregularities  (e.g.,  screen  support  spokes).  Figure  17  shows 
typical  data  from  the  pressure  transducer  nearest  the  leading  edge  on 
the  suction  surface  after  each  stage  of  processing.  The  top  curve  is 
the  raw  data  over  one  cycle  after  scaling  to  proper  coefficient  form. 
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The  middle  curve  is  the  result  of  averaging  the  signal  over  five 
cycles,  and  the  bottom  curve  is  the  same  data  after  applying  the 
digital  filtering. 

The  pressure  coefficients  are  integrated  to  obtain  normal  force 
and  pitching  moment  coefficients  (Ref.  21).  The  sign  conventions 
used  ih  presenting  these  parameters  are  illustrated  in  Fig.  18.  The 
normal  force  coefficient  is  given  by 


V "i 


(13) 


where  the  w^'s  are  the  Gaussian  weighting  coefficients.  Transducer 
numbers  1-6  are  located  on  the  suction  surface  and  transducer  numbers 
7-12  are  located  on  the  pressure  surface  at  the  same  respective 


FIG.  18  SIGN  CONVENTIONS  FOR  NORMAL  FORCE  AND 
PITCHING  MOMENT  COEFFICIENTS 
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chordwise  positions.  The  pitching  moment  coefficient  about  the 
midchord  is  given  by 


C = I (C  ) w. (x.  - 0.5) 

M j p . . 1 1 


(14) 


i=l  Fi+6 

where  the  x^'s  are  the  chordwise  positions  of  the  measuring  stations 
expressed  as  fractions  of  the  chord.  The  Gaussian  positions  and 
weighting  factors  used  in  these  equations  are  presented  in  Table  I. 
Even  with  only  six  chordwise  stations,  the  Gaussian  integration 
technique  may  be  expected  to  yield  normal  force  and  pitching  moment 
coefficients  accurate  to  within  three  percent  (Ref.  22). 


Transient  Surface  Flow  Phenomena 
Hot  film  skin  friction  gages  were  positioned  adjacent  to  each 
pressure  transducer  in  an  attempt  to  correlate  the  pressure  fluctua- 
tions with  surface  flow  phenomena.  The  hot  film  gage  is  a resistance 
element  which  is  heated  by  passing  electrical  current  through  it.  A 
constant-temperature  anemometer  circuit  is  employed  to  keep  its 
resistance  and  hence  its  thermometer  constant.  By  monitoring  the 
power  necessary  to  accomplish  this,  it  is  possible  to  obtain  a 
qualitative  measure  of  the  heat  transfer  to  the  fluid.  McCroskey 
(Ref.  28)  gives  this  relation  as 

I2R  a A + Q , (15) 

where  I is  the  current  through  the  gage  and  R is  its  resistance.  The 


heat  transfer  is  composed  of  two  parts;  the  nonflow  heat  losses,  A, 
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and  a nondimens ional  forced  convection  loss,  Q.  Normally,  the 
current  through  the  gage  is  not  measured  directly.  Instead,  the 
anemometer  output  voltage  is  measured;  however,  the  relationship 
between  the  gage  current  and  output  voltage  is  such  that 

E2  oc  A + Q (16) 

This  voltage  is  composed  of  a steady  part,  E,  and  a time-varying 
part,  e,  which  is  much  less  than  E.  For  the  present  study,  only  quali- 
tative data  are  desired  from  the  skin  friction  gages;  hence,  only  the 
time  varying  part  of  the  heat  transfer  need  be  retained.  Neglecting 
higher  order  terms  in  the  time-varying  voltage  yields 

e a Q 07) 

/ 

where  Q is  the  time  varying  portion  of  the  forced  convection.  From 
Ref.  29 

Q a T . (18) 

w 

where  tw  is  the  shearing  stress  at  the  wall.  The  skin  friction 
coefficient  is  defined  as  (Ref.  29) 


(19) 


Considering  only  the  time-varying  part  of  the  skin  friction  coefficient, 
we  may  write 


- 4 


e (— ) 
lWi 


(20) 
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I 
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Although  the  quantity  on  the  right  of  Eq.  (20)  has  the  units  of 
volts,  it  has  been  corrected  for  local  velocity  variations  and  will 
be  used  as  an  indicator  of  the  variations  in  the  skin  friction 
coefficient.  It  is  not  valid  to  compare  the  amplitudes  of  this 
parameter  at  different  flow  conditions  because  the  constant  of 
proportionality  need  not  be  the  same  in  the  two  cases.  Furthermore, 

only  the  varying  part  of  the  skin  friction  is  being  considered;  thus, 

^1  2/3 

when  the  quantity  eC^-ji)  ' becomes  zero,  this  does  not  mean  that 

W1 

the  total  skin  friction  coefficient,  C^,  is  zero.  In  presenting 
this  data,  no  averaging  or  filtering  techniques  are  employed  because 
the  high  frequency  fluctuations  in  the  parameter  defined  by  Eq.  (20) 
are  needed  to  judge  whether  transition  or  separation  has  occurred. 


Velocity  Measurements  at  the  Rotor  Exit  Plane 
A complete  description  of  the  calibration  procedure  developed 
for  the  triaxial  hot-film  probe  used  in  this  experiment  appears  in 
Appendix  V.  It  was  determined  that  the  response  of  each  of  the 
sensors  to  a purely  normal  velocity  obeyed  the  relation 

E2  - E2  - k^V  (Ts  - T„)  (21) 

if  one  postulated  that  the  empirical  constants  Eq  and  k might  have 
one  value  at  high  velocities  and  another  at  low  velocities.  Here,  E 
is  the  bridge  output  voltage  and  p and  V are  the  fluid  density  and 
velocity  (assumed  to  be  normal  to  the  sensor).  The  temperatures  Tg 


and  Tb  are  those  of  the  sensor  and  free  stream,  respectively.  The 


51 


constants  Eq  and  k were  determined  in  the  calibration  procedure 
(Appendix  V).  The  equation  may  be  rewritten  to  give  an  effective 
normal  velocity  as  a function  of  these  calibration  constants,  the 
fluid  properties,  and  the  anemometer  output  voltage. 


q - [- 


E2  - Eo2 


(Ts  - TJ 


-] 


(22) 


It  was  found  that  the  best  relation  between  the  sensed  normal  velocity, 
q,  and  the  actual  normal  velocity,  Un,  was  given  by 


q2  * U2  f (0,  4>)  . 


(23) 


The  function  f(9,  ♦)  is  dependent  on  the  pitch  and  yaw  angles  of  the 
flow  relative  to  the  probe  (9  and  *,  respectively)  and  is  evaluated 
using  the  calibration  data.  This  function  is  different  for  each 
sensor.  Given  values  of  0 and  ♦,  one  may  evaluate. the  function  for 
each  wire.  Making  use  of  the  fact  that  the  three  sensors  are  mutually 
orthogonal,  one  may  write  Eq.  (23)  for  each  sensor. 

2 


V-  = U ♦ U 
t 1 2 


(24) 


n 


h 


U - + u 
1 3 


u 2 + u 2 


(25) 


(26) 
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In  these  equations,  is  the  value  of  the  function  f(e,  ♦)  for  the 
^th  sensor  evaluated  at  the  assumed  angles  6 and  ♦ and  is  the 

velocity  parallel  to  the  i**1  sensor.  These  equations  may  be  solved 
easily  for  the  squared  velocity  components  u^.  Since  there  is 
no  way  to  evaluate  the  algebraic  sign  of  u^  form  this  process,  it 
must  be  known  from  some  other  source.  In  the  present  experiment,  the 
orientation  of  the  probe  was  such  that  the  velocity  components  may  be 
assumed  to  have  the  positive  sign.  Using  these  velocities,  the 
assumed  angles  were  updated  and  the  process  repeated  until  convergence 
was  obtained.  A more  complete  description  of  the  iteration  process 
and  the  equations  involved  may  be  found  in  Appendix  V.  The  final 
result  is  expressed  in  terms  of  the  relative  velocity  at  the  exit 
plane,  , the  flow  angle  in  the  axial-circumferential  plane,  Bj, 
and  the  flow  angle  in  the  axial-radial  plane,  If  one  defines  the 

velocities  Vj  , V2,  v-j  to  be  the  radial,  circumferential,  and 
axial  velocity  components,  respectively,  then  relative  velocity  and 
flow  angles  are  given  by: 
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The  relative  velocity  and  the  two  flow  angles  were  computed  for 
five  revolutions  with  roughly  200  data  points  per  cycle.  These  data 
were  represented  by  Fourier  Series  in  terms  of  the  circumferential 
position  in  the  flow  field  so  that  the  parameters  of  interest  could 


later  be  reconstructed  at  any  desired  rotor  position. 
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RESULTS 

The  notation  used  herein  is  that  employed  in  turbomachine 
work  by  UTRC.  In  particular,  the  variable  against  which  most  of 
the  data  are  plotted  is  the  rotor  inlet  angle,  8^.  This  angle  has 
previously  been  defined  in  Eq.  4 (see  Fig.  6).  It  should  be  noted 
that  increasing  Bj  corresponds  to  decreasing  angle  of  attack. 

Further,  the  mechanism  by  which  the  inlet  angle  is  changed  is  to  vary 
the  axial  velocity  at  constant  rotational  speed.  Consequently,  the 
effects  of  the  varying  inlet  angle  (or  angle  of  attack)  cannot  be 
totally  separated  from  the  effects  of  the  varying  velocity.  Proper 
nondimensionalizat ion  of  the  various  measured  quanities  can  eliminate 
the  velocity  effects  in  the  case  of  quasi-steady  flow.  However, 
unsteady  effects  of  the  varying  velocity  are  still  present.  The 
method  employed  to  illustrate  the  effects  of  unsteadiness  in  the  flow 
is  to  plot  the  unsteady  response  of  the  rotor  undergoing  inlet 
distortion  superimposed  on  its  steady-state  behavior.  This  steady- 
state  behavior  was  determined  by  running  the  rotor  at  different  inlet 
angles  in  undistorted  flow. 

Overall  Performance 

The  annulus  averaged  pressure  rise  coefficient  across  the  rotor 
(defined  by  Eq.  (5))  is  shown  in  Fig.  19  for  both  undiscorted  and  dis- 
torted flow.  For  distorted  flow,  the  pressure  rise  is  circumferentially 
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averaged  and  plotted  as  a function  of  the  average  inlet  angle.  The 
upper  branch  of  each  curve  represents  unstalled  flow  and  the  lower 
branch  represents  stalled  operation.  (Refer  back  to  Fig.  3 and  the 
associated  text  for  a more  complete  discussion  of  the  discontinuous 
behavior  of  the  curves.)  The  undistorted  curve  shows  a significantly 
greater  peak  pressure  rise  and  is  slightly  higher  than  the  distorted 
curve  over  the  higher  inlet  angles.  However,  at  the  lower  inlet 
angles,  where  the  blade  loading  is  higher,  the  unsteady  effects  cause 
the  distorted  flow  condition  to  give  slightly  higher  pressure  rise. 
Note  that  the  incidence  angle  (relative  to  the  mean  camber  line)  at 

which  stall  occurs  is  very  low,  only  1.5°  for  undistorted  flow.  (The 

• • • • , 
incidence  angle  is  the  difference  between  the  stagger  angle,  g[ , which 

is  30°  in  this  case  and  the  inlet  angle.)  This  is  possible  because 

the  blade  profile  is  highly  cambered  and  can  produce  large  turning 

angles  (and  consequent  high  normal  forces)  at  low  incidence  angles. 

The  low  incidence  angle  at  stall  is  common  to  highly  cambered  blades. 

The  loss  coefficient  of  Eq . (6)  is  plotted  as  a function  of 

inlet  angle  in  Fig.  20  for  undistorted  flow  only.  Suitable  equipment 

for  measuring  the  unsteady  loss  was  not  available  when  this  experiment 

was  run. 

To  complete  the  presentation  of  the  overall  performance  para- 
meters, we  will  consider  the  variation  of  the  exit  flow  angle  as  a 
function  of  the  inlet  angle  (again,  for  undistorted  flow  only)  as 
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shown  in  Fig.  21.  Curves  are  presented  for  the  basic  wedge  probe 
measurements,  for  a momentum  theory  calculation  (see  Appendix  V),  and 
for  cascade  data  taken  from  Ref.  24  for  a similar  profile  having  the 
same  stagger  angle  and  gap-to-chord  ratio.  The  three  curves  show 
reasonable  agreement  for  only  a narrow  range  of  inlet  angles  (39°  < 

8^  < 45*)  with  large  differences  at  the  lower  values  of  inlet  angle 
(high  blade  loading).  Several  factors  may  account  for  these  differ- 
ences and  most  of  them  are  in  some  way  associated  with  the  three- 
dimensional  lity  of  the  flow.  First,  the  correction  for  axial  velocity 
variation  across  the  rotor  i6  empirical  and  may  not  be  accurate  at 
high  blade  loadings.  This  would  introduce  errors  into  the  exit  angle 
as  calculated  from  the  wedge  probe  measurements.  Second,  the  radial 
velocity  may  not  be  zero  as  was  assumed.  The  midspan  behavior  would 
not  be  representative  of  the  average  blade  performance;  hence,  the 
wedge  probe  measurements  could  not  be  expected  to  agree  with  the 
two-dimensional  cascade  data.  Furthermore,  the  radial  equilibrium 
correction  to  the  rotor  pressure  rise  would  be  in  error  thus  intro- 
ducing error  into  the  momentum  theory  calculations.  Finally,  the 
blade  wakes  may  introduce  enough  unsteadiness  into  the  flow  so  that 
the  flow  angles  measured  by  the  wedge  probes  are  in  error  (c.f.,  Ref. 
18).  At  the  higher  inlet  angles,  the  exit  angle  as  calculated  from 
the  wedge  probe  measurements  is  in  good  agreement  with  the  cascade 
data  but  that  from  momentum  theory  appears  to  be  diverging.  Since 
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the  momentum  theory  calculations  is  subject  to  the  greatest  number  of 
errors,  the  observed  differences  are  not  considered  excessive.  It 
must  be  surmised  that  the  placing  of  a single  wedge  probe  at  the 
midspan  location  is  not  an  adequate  means  of  measuring  the  exit 
angle.  A more  prudent  method  would  be  to  use  a spanwise  array  of 
high  response  (hot-wire  or  hot-film)  probes. 

Analysis  of  Unsteady  Data 

It  has  been  stated  that  the  method  used  to  determine  the  effects 
of  the  unsteadiness  generated  by  the  distortion  will  be  to  plot  the 
various  parameters  as  a function  of  the  inlet  angle.  This  process 
normally  results  in  a loop  which  will  generally  be  superimposed  on 
the  quasi-steady  characteristic.  Before  discussing  the  experimental 
data,  several  simple  examples  will  be  given  to  illustrate  how  these 
loops  are  derived  and  how  their  character  may  be  related  to  changes 
in  the  amplitude  and  phase  angle  of  the  response.  Consider  the  two 
variables  x and  y which  are  sinusoidal  functions  of  time. 

x = A sin  <ot  (30) 

y = B sin  (to  + <J>)  (31) 

Suppose  x is  a displacement  and  y is  a response  due  to  the  displace- 
ment. The  amplitude  of  the  displacement  is  A,  that  of  the  response 
is  B,  and  $ is  the  phase  angle  by  which  the  response  leads  the  dis- 
placement. Some  ways  in  which  x and  y may  vary  are  shown  in  Fig.  22 
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together  with  the  loops  which  result  from  plotting  y as  a function  of 
x.  If  the  two  variables  are  in  phase,  y is  a linear  function  of  x 
and  the  slope  is  indicative  of  the  relative  amplitudes  of  x and  y 
(e.g.,  solid  and  dashed  lines  in  Fig.  22).  If  the  two  variables  are 
not  in  phase  plotting  y as  a function  of  x produces  an  ellipse 
(dash-dot  lines  in  Fig.  22).  The  small  arrows  indicate  the  direction 
in  which  the  loop  is  traversed.  When  the  slope  of  the  major  axis  of 
the  loop  is  positive,  the  counterclockwise  sense  shown  in  Fig.  22 
indicates  that  the  response  lags  the  displacement.  If  the  slope  of 
the  loop  is  negative,  a clockwise  sense  indicates  a lagging  response 
and  a counterclockwise  sense  indicates  a leading  response. 

In  the  case  of  a compressor  rotor,  the  inlet  angle  may  be 
considered  to  be  the  displacement  mentioned  above  and  the  blade 
pressures,  etc.  are  the  responses.  Acutally,  none  of  these  are  pure 
sinusoids  and  the  generated  loops  will  not  be  simple  ellipses. 

However,  certain  basic  similarities  do  exist  between  this  illustrative 
example  and  the  actual  measured  data.  Therefore,  comparing  the 
unsteady  responses  (obtained  while  operating  under  the  influence  of 
the  inlet  distortion)  to  their  steady  counterparts  (obtained  with 
undistorted  flow)  will  reveal  the  nature  and  extent  of  unsteady 
effects.  First,  the  slope  of  the  loop  relative  to  the  steady  state 
curve  will  indicate  whether  unsteadiness  has  changed  the  amplitude  of 
the  response.  Second,  the  openness  of  the  loop  will  show  phase  shift 
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and,  if  the  loop  has  a definite  positive  or  negative  slope,  the 
direction  in  which  the  loop  is  traversed  will  indicate  whether  the 
response  leads  or  lags  the  inlet  angle.  A loop  which  is  horizontal 
indicates  a phase  angle  near  90*  (or  270*).  It  should  be  noted  that 
in  the  inlet  distortion  problem  the  leading  and  lagging  phase  angles 
do  not  imply  stability  or  instability  as  is  the  case  with  an  oscil- 
lating airfoil  (or  a cascade  of  airfoils);  rather,  the  phase  angle  is 
simply  an  indication  of  unsteadiness. 

Blade  Pressures 

The  most  instructive  form  for  presenting  the  blade  pressures 
is  the  difference  between  the  pressure  coefficient  on  the  pressure 
surface  and  that  on  the  suction  surface.  This  differential  pressure 
coefficient  is  given  by 


\ ’ [‘ 


- |Cp  - Cp 

press  suet 


L 


(32) 


where  the  subscript  i refers  to  the  chordwise  station  in  question. 
Since  the  individual  coefficients  were  already  represented  as  Fourier 
Series  (of  the  form  given  by  Eq.  (12)),  it  is  a simple  matter  to  ob- 
tain the  series  representations  of  the  differential  pressure  coeffi- 
cients. These  were  converted  to  an  amplitude/phase  angle  format 


AC_  ■ a + I a cos  (nGH-$  ) 
o n n 


(33) 


where  aQ  is  the  steady  component  and  an  and  +n  are  the 
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amplitude  and  phase  angle  of  the  nth  harmonic,  respectively.  The 
upper  limit  of  the  summation,  N,  is  a function  of  the  cutoff  frequency 
of  the  digital  filter  and  8 is  the  angular  position  in  the  distortion 
field  with  the  point  of  maximum  axial  velocity  being  8 “ 0. 

Three  flow  conditions  were  chosen  as  representative  of  those 
for  which  data  were  available.  The  mean  inlet  angles  associated  with 
these  flow  conditions  are  indicated  by  tic  marks  on  the  undistorted 
rotor  pressure  rise  curve  (Fig.  23).  These  angles  represent  low 
(7.15),  medium  (7.9),  and  high  (7.4)  blade  loading.  It  is  instructive 
to  compare  the  amplitudes  and  phase  angles  of  the  first  few  harmonics 
of  the  differential  pressure  coefficients  and  the  integrated  normal 
force  and  pitching  moment  coefficients  (Table  II).  One  concludes 
that  the  first  harmonic  is  dominant  (as  much  as  an  order  of  magnitude 
greater  than  the  higher  harmonics).  The  response  at  the  leading  edge 
stations  is  significantly  greater  than  that  at  stations  farther  aft 
on  the  blade.  Tms  peaking  of  the  response  at  the  leading  edge  seems 
to  be  more  pronounced  at  higher  blade  loading.  The  lack  of  activity 
near  the  trailing  edge  is  also  noteworthy. 

The  differential  pressure  coefficient  loops  for  each  of  the  six 
chordwise  measuring  stations  are  presented  in  Fig.  24.  At  the 
leading  edge  (x/c  “ 0.034),  one  observes  that  the  amplitude  of  the 
response  has  been  increased  and  that  the  response  consistently  lags 

4 

the  inlet  angle.  Horlock,  Greitzer,  and  Henderson  (Ref.  30)  have 

D 

(I 


TABLE  II 
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HARMONIC  CONTENT  OF  DIFFERENTIAL  PRESSURE, 
NORMAL  FORCE  AND  P TCHING  MOMENT  COEFFICIENTS. 


(a)  Point  7.4  Bj  = 28.9  + .061  cos  6 
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CM 

0.052 

178.5° 

0.005 

177.5 

a3 

♦3 

a4 

*4 

x/c  = 0.034 

0.129 

289.1° 

0.042 

283.0 

x/c  = 0.169 

0.007 

214.8° 

0.002 

52.8 

x/c  = 0.381 

0.018 

65.2° 

0.007 

22.1' 

x/c  = 0.619 

0.019 

251.3° 

0.006 

58.2' 

x/c  = 0.831 

0.006 

91.1° 

0.007 

309.8' 

x/c  = 0.966 

0.007 

331.7° 

0.006 

309.2' 

CN 

0.011 

288.1° 

0.005 

328.7' 

CM 

0.005 

288.9° 

0.001 

276.4' 

TABLE  II  (Cont'd) 
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HARMONIC  CONTENT  OF  DIFFERENTIAL  PRESSURE, 
NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS. 


(b)  Point. 7. 9 Bj  = 36.9  + .068  cos  6 


al 

n 

a2 

♦2 

x/c  = 0.034 

0.737 

173.9° 

0.009 

296.3 

x/c  = 0.169 

0.244 

181.2° 

0.027 

204.3 

x/c  = 0.381 

0.074 

252.0° 

0.029 

271  .7 

x/c  = 0.619 

0.206 

317.4° 

0.039 

288.5 

x/c  = 0.831 

0.030 

271.9° 

0.036 

131.6 

x/c  = 0.966 

0.133 

173.7° 

0.023 

109.8 

CN 

0.101 

207.8° 

0.013 

242.4 

°M 

0.045 

172.2° 

0.003 

271  .0 

a3 

$3 

a4 

♦4 

x/c  = 0.034 

0.063 

166.0° 

0.039 

114.9 

x/c  = 0.169 

0.037 

136.6° 

0.008 

109.8 

x/c  = 0.381 

0.002 

305.2° 

0.004 

75.6 

x/c  = 0.619 

0.020 

19.6° 

0.017 

33.2 

x/c  = 0.831 

0.012 

53.9° 

0.004 

115.7' 

x/c  = 0.966 

0.016 

115.3° 

0.011 

190.7' 

CN 

0.01 1 

115.0° 

0.008 

86.7' 

CM 

0.005 

171.1° 

0.002 

112.8' 

TABLE  II  (Concluded) 
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HARMONIC  CONTENT  OF  DIFFERENTIAL  PRESSURE, 
NORMAL  FORCE  AND  PITCHING  MOMENT  COEFFICIENTS. 


(c)  Point. 7. 15  = 42.4  + .064  cos  e 


al 

♦l 

a2 

*2 

x/c  = 0.034 

0.858 

173.0° 

0.053 

265.0° 

x/ c = 0.169 

0.312 

172  .4° 

0.028 

285.8° 

x/c  = 0.381 
x/c  = 0.61S 

0.092 

188.1° 

312.8° 

0.013 

248.5° 

x/c  = 0.831 

0.064 

321.0° 

0.006 

222.7° 

x/c  = 0.966 

0.102 

206.6° 

0.019 

117.4° 

CN 

0.136 

185.3° 

0.015 

272  .9° 

°M 

0.057 

167.7° 

0.004 

276.4° 

a3 

♦3 

a4 

♦4 

x/c  = 0.034 

0.140 

150.3° 

0.013 

125.1° 

x/c  = 0.169 

0.021 

158.4° 

0.008 

357.4° 

x/c  = 0.381 
x/c  = 0.619 

0.008 

190.0° 

0.007 

213.1° 

x/c  = 0.831 

0.007 

250.9° 

0.007 

216.7° 

x/c  = 0.966 

0.011 

187.3° 

0.009 

60.2° 

Cn 

0.017 

163.2° 

0.004 

248.0° 

CM 

0.007 

151  .7° 

0.001 

71.5° 
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studied  the  response  of  a cascade  to  low  frequency  inlet  distortions 
analytically.  Their  results  show  that  the  magnitude  of  the  unsteady 
lift  coefficient  increases  with  increasing  reduced  frequency  and  that 
the  amplification  effect  exists  even  at  the  low  reduced  frequencies 
(k  = 0.1)  encountered  in  the  present  experimental  program.  Although 
Horlock,  et  al . , do  not  calculate  the  behavior  of  the  surface  pressures 
explicitly,  it  is  reasonable  to  infer  that  their  variations  would 
also  be  amplified.  Thus,  the  increased  amplitude  in  the  leading  edge 
pressure  coefficients  is  not  totally  unprecedented.  The  blade 
penetrates  considerably  beyond  the  steady-state  stall  point,  the  left 
end  of  the  dashed  line,  resulting  in  a pressure  differential  far  in 
excess  of  that  possible  in  undistorted  flow.  The  phenomenon  of  stall 
delay  on  an  isolated  airfoil  oscillating  in  pitch  is  widely  recognized 
(Refs.  31-33).  This  delay  is  a strong  function  of  reduced  frequency 
and  is  normally  not  so  pronounced  at  the  low  reduced  frequencies 
encountered  in  this  experiment.  Further,  this  delay  is  commonly 
attributed  to  a decrease  in  the  amplitude  of  the  normal  force  variation 
as  reduced  frequency  increases  (Refs.  32,  34,  and  35).  The  fact  that 
the  amplitude  of  the  response  is  increased  in  the  present  experiment 
and  in  the  theoretical  analysis  of  Ref.  30  must  be  attributed  to 
multiblade  effects.  One  manner  in  which  multiblade  effects  may 
manifest  themselves  is  in  the  existence  of  an  interblade  phase  angle. 

In  the  present  experiment,  the  interblade  phase  angle,  o,  was  - 12.9*. 
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(This  is  obtained  from  a stationary  360°  wave  divided  over  28  blades 
so  that  it  appears  to  be  a backward-travelling  wave.)  An  experimental 
investigation  in  a cascade  of  airfoils  oscillating  in  pitch  (Ref.  36) 
has  shown  that  the  interblade  phase  angle  has  a strong  influence  on 
the  amplitudes  and  phase  angles  of  the  normal  force  and  pitching 
moment  coefficients.  It  is  thus  reasonable  to  suppose  that  the 
interblade  phase  angle  present  in  the  inlet  distortion  experiment 
is  at  least  partially  responsible  for  the  marked  departure  from  the 
steady  state  behavior  as  evidenced  by  the  change  in  the  amplitude  of 
the  response  and  the  pronounced  stall  delay. 

The  pressure  loops  exhibit  fluctuations  which  are  greater  at  the 
lower  values  of  for  each  loop.  These  fluctuations  occur  at 
exactly  half  the  screen-segment  passing  frequency.  The  screen  which 
produces  the  distortion  is  constructed  of  segments  of  wire  mesh  of 
different  porosities.  Therefore,  the  axial  velocity  variation  is  not 
perfectly  smooth.  This  may  be  thought  of  as  a higher  order  velocity 
variation  superimposed  on  the  once-per-revolution  variation.  It  thus 
appears  that  this  higher  frequency  disturbance  is  inducing  a (rela- 
tively) larger  response.  During  the  experiment,  an  overload  occurred 
in  the  static  data  system  requiring  that  the  average  pressure  mea- 
surements for  x/c  “ 0.034  be  corrected  as  described  in  Appendix 


VII.  No  other  stations  were  affected. 


I 

I 

I 
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At  x/c  = 0.169  the  response  is  much  less  than  at  x/c  = 0.034, 
although  it  is  still  slightly  amplified  over  the  undistorted  curve. 

For  the  high  blade  loading  of  point  7.4,  the  response  leads  the 
input.  For  point  7.9,  the  response  leads  at  the  higher  instantaneous 
inlet  angles  and  lags  at  the  lower  angles.  For  the  lowest  blade 
loading,  point  7.15,  the  response  lags  the  inlet  angle  over  most  of 
the  cycle.  The  phase  angle  of  the  pressure  response  relative  to  the 
inlet  angle  thus  decreases  as  the  inlet  angle  increases. 

At  x/c  = 0.381  the  loop  representing  point  7.4  is  more  open 
than  the  other  two.  This  is  because  the  first  harmonic  amplitude  is 
greater  for  this  point  and  the  phase  angle  is  quite  near  270°. 

(Phase  angles  of  90  or  270  degrees  will  produce  the  most  open  loops.) 
As  with  x/c  = 0.169,  the  phase  angle  consistently  decreases  as  the 
average  inlet  angle  increases.  The  loops  have  all  become  considerably 
smoother,  perhaps  indicating  that  the  fluctuations  appearing  at  the 
leading  edge  have  been  damped  out  by  the  presence  of  the  neighboring 
blades . 

All  the  loops  for  x/c  = 0.619  are  counterclock-wise,  with  the 
loop  for  point  7.9  being  more  open.  The  phase  angle  of  the  pressure 
response  does  not  decrease  with  increasing  inlet  angle  as  was  previ- 
ously observed;  instead,  it  remains  relatively  constant.  The  pressure 
fluctuations  first  observed  at  the  leading  edge  have  been  further 
reduced  in  amplitude.  The  loops  for  x/c  = 0.831  and  x/c  “ 0.966  show 
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very  little  activity.  The  undistorted  characteristic  is  very  nearly 
horizontal  in  both  cases  indicating  that  the  pressure  coefficients  at 
these  stations  are  essentially  independent  of  inlet  angle.  The 
pressure  coefficient  difference  appears  to  be  approaching  zero  at  the 
trailing  edge . 

The  variations  in  the  amplitude  of  these  responses  were  well 
illustrated  by  Fig.  24;  however,  the  phase  angle  variations  are  not 
so  obvious.  The  phase  angle  is  a strong  indicator  of  unsteadiness, 
perhaps  more  so  than  is  amplitude.  (See,  for  example,  Ref.  34  or 
Ref.  33.)  The  phase  angle  of  the  differential  pressure  coefficients 
and  the  integrated  normal  force  and  pitching  moment  coefficients  are 
presented  in  Table  III.  The  phase  angles  for  the  two  forward  measur- 
ing stations  and  for  the  normal  force  and  pitching  moment  are  shown 
as  a function  of  the  mean  inlet  angle  in  Fig.  25.  Because  the 
reference  is  to  the  inlet  angle,  8j,  rather  than  to  an  incidence 
angle,  quasi-steady  behavior  would  be  indicated  by  a phase  angle  of 
180'  rather  than  0°  as  is  normally  the  case.  The  phase  angle  is  seen 
to  vary  smoothly  with  the  mean  inlet  angle  indicating  that  the 
unsteadiness  of  the  response  is  a strong  function  of  blade  loading. 

To  complete  the  discussion  of  blade  pressures,  we  will  now 
consider  the  pressure  distribution  across  the  chord  of  the  blades. 
This  is  simply  another  method  of  presenting  the  pressure  data  to 
emphasize  the  way  in  which  the  pressures  at  different  positions  on 
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TABLE  III 

FIRST  HARMONIC  PHASE  ANGLES  OF  THE 
DIFFERENTIAL  PRESSURES,  NORMAL  FORCE,  AND  PITCHING  MOMENT 


Chordwi se 

St  at  ion 

28.9 

33.5 

1 

36.9 

40.4 

42.2 

45.2 

0.034 

172* 

177* 

174* 

175* 

173* 

172* 

0.169 

197* 

188* 

181* 

178* 

172* 

169* 

0.381 

285* 

282* 

252* 

219* 

188* 

173* 

0.619 

262* 

319* 

317* 

316* 

313* 

319° 

0.831 

164° 

162* 

272* 

306* 

321* 

330* 

0.966 

111' 

172° 

174* 

191* 

207° 

220* 

CN 

201* 

199° 

208* 

204* 

185* 

168* 

179* 

185* 

172* 

169* 

168* 

166* 

FIG.  25  PRESSURE,  FORCE,  AND  MOMENT  PHASE  ANGLE  AS  A FUNCTION 
OF  AVERAGE  INLET  ANGLE 
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the  blade  relate  to  one  another  at  a given  inlet  angle.  Four  inlet 
angles  will  be  considered  for  each  flow  condition.  These  will  be  the 
minimum  inlet  angle  ( 9 = 180*),  the  mean  inlet  angle  with  the  inlet 
angle  increasing  ( 9 * 270°),  the  maximum  inlet  angle  (9  = 0*),  and 
the  mean  inlet  angle  with  the  inlet  angle  decreasing.  The  pressure 
coefficient  distributions  for  the  suction  surface  (denoted  by  an  S) 
and  the  pressure  surface  (denoted  by  a P)  are  presented  for  the  three 
representative  flow  conditions  in  Figs.  26  through  28,  respectively. 
Although  the  distribution  of  measuring  stations  gives  adequate  accur- 
acy in  the  integration  to  obtain  normal  force  and  pitching  moment 
coefficients,  it  is  not  sufficiently  dense  near  the  leading  edge  to 
define  the  shape  of  the  pressure  profile  at  the  leading  edge.  In 
particular,  it  is  not  possible  to  resolve  the  location  or  magnitude 
of  the  leading  edge  suction  peak  under  high  load. 

It  must  be  remembered  that  decreasing  inlet  angle  corresponds 
to  increasing  blade  incidence.  Figure  26  for  point  7.4  shows  a 
relatively  sharp  peaking  of  the  leading  edge  suction  (although  the 
exact  value  of  the  peak  and  its  precise  location  cannot  be  determined) 
which  is  most  pronounced  at  the  minimum  instantaneous  inlet  angle 
(maximum  incidence).  The  peak  is  in  evidence,  although  not  as 
strong,  for  all  inlet  angles  at  this  flow  condition.  As  the  mean 
inlet  angle  is  increased  to  the  flow  condition  of  point  7.9  (Fig. 

27),  this  peaking  is  lessened.  In  fact,  the  peak  at  the  extreme 


(P  DENOTES  PRESSURE  SURFACE.  S DENOTES  SUCTION  SURFACE) 


FIG.  26  PRESSURE  COEFFICIENT  DISTRIBUTIONS  FOR  POINT  7.4 


leading  edge  is  not  present  at  the  maximum  instantaneous  inlet  angle 
for  this  flow  condition.  Rather,  the  suction  surface  pressure 
distribution  for  this  inlet  angle  is  smoother  with  a mild  peaking 
farther  aft  on  the  blade.  This  trend  continues  with  an  additional 
increase  in  mean  incidence  angle  to  attain  the  flow  condition  of 
point  7.15  (Fig.  28).  The  leading  edge  suction  peak  is  evident  only 
for  the  lowest  instantaneous  inlet  angle.  Furthermore,  at  the 
maximum  inlet  angle  for  this  flow  condition,  the  pressure  on  the 
suction  surface  is  actually  higher  than  that  on  the  pressure  surface 
This  pressure  reversal  is  due  to  the  extreme  negative  incidence 
(omc^  = -15.4°)  at  this  inlet  angle.  It  may  also  be  noted  that  the 
variations  in  the  pressure  distributions  occur  primarily  in  the 
leading  edge  region  over  the  entire  range  of  inlet  angles  in  these 
three  figures.  Such  behavior  is  fairly  typical  of  cascaded  airfoils 
Normal  Force  and  Pitching  Moment  Coefficients 

The  normal  force  and  pitching  moment  coefficients  were  obtained 
by  integrating  the  pressure  coefficient  distributions.  It  is  thus 
appropriate  to  relate  the  behavior  of  normal  force  and  pitching 
moment  coefficient  loops  to  the  pressure  coefficient  loops.  The 
normal  force  loop  (Fig.  29)  is  dominated  by  the  pressures  at  x/c  = 
0.381  and  x/c  = 0.619  because  these  two  stations  have  the  highest 
weighting  coefficients  in  the  Gaussian  integration  scheme.  (See  Eq. 
9 and  Table  I).  Since  all  three  pressure  loops  at  both  of  these 
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FIG.  29  NORMAL  FORCE  COEFFICIENT  AS  A FUNCTION  ‘ 
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stations  are  counterclockwise  (see  Fig.  24),  it  is  not  surprising 
that  the  normal  force  loops  are  also  counterclockwise.  Even  the  fact 
that  the  loop  for  point  7.5  is  almost  completely  collapsed  on  itself 
can  be  traced  directly  to  these  pressure  loops.  At  first,  it  appears 
that  the  normal  force  loops  lie  further  above  the  quasi-steady  curve 
than  would  be  expected  from  the  pressure  loops.  This  illusion  is  due 
to  the  expanded  scale  used  to  plot  the  normal  force  loops.  The  large 
fluctuations  observed  in  the  pressure  loops  for  the  leading  edge  sta- 
tions are  present  and  have  also  been  magnified  by  the  scale  change. 
Perhaps  the  most  interesting  point  is  that  the  maximum  normal  force 
coefficient  in  the  distorted  flow  is  roughly  25%  greater  than  the 
maximum  quasi-steady  value  and  is  due  in  part  to  the  excursion  past 
the  steady-state  stall  point.  Of  course,  the  variation  in  actual 
load  is  not  nearly  so  great  since  this  high  normal  force  coefficient 
occurs  in  the  low  velocity  region  of  the  distortion. 

The  pitching  moment  loop  (Fig.  30)  is  virtually  identical 
to  the  pressure  loop  at  x/c  = 0.034.  Although  the  stations  at  x/c  * 
0.0381  and  x/c  = 0.619  have  the  highest  weighting  coefficients,  they 
are  located  symmetrically  about  the  moment  center  (mid-chord)  and 
effectively  cancel  each  other.  (See  Eq.  10  and  Table  I.)  The 
pitching  moment  is  thus  dominated  by  the  strong  leading  edge  response. 
The  loops  are  all  basically  clockwise,  indicating  that  the  pitching 
moment  lags  the  inlet  angle,  and  becomes  more  open  as  the  average 
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FIG.  30  PITCHING  MOMENT  COEFFICIENT  AS  A FUNCTION  OF  INLET  ANGLE 


r 


r 


I 

I 

I 

! 

I 


I 

I 

I 

I 

I 


8A 

inlet  angle  increases,  indicating  a greater  lag  for  the  lighter  blade 
loading.  This  behavior  was  also  evident  in  Fig.  25. 

Comparison  to  Cascade  Data 

The  aerodynamics  of  a rotor  undergoing  inlet  distortion  and  a 
cascade  oscillating  in  pitch  are  not  identical.  However,  if  one 
compares  the  phase  plane  diagrams  of  the  Sears  function  for  a convec- 
ted  sinusoidal  gust  (Ref.  34)  and  the  Theodorsen  function  for  a 
sinusoidally  oscillating  airfoil  (Ref.  35),  it  is  apparent  that  the 
convected  gust  problem  and  the  oscillating  airfoil  problem  do  have 
similar  solutions  at  low  reduced  frequencies.  It  is  thus  logical  to 
expect  the  cascade  tests  of  Ref.  36  to  reproduce  the  basic  trends 
observed  in  the  inlet  distortion  experiment.  This  is  indeed  the 

case.  The  peaking  of  the  suction  near  the  leading  edge  and  the 

/ 

concentration  of  the  pressure  variations  near  the  leading  edge 
observed  in  the  present  experiment  are  duplicated  in  the  cascade 
experiment.  At  an  equivalent  reduced  frequency  and  zero  interblade 
phase  angle,  the  normal  force  coefficient  for  the  cascade  generally 
leads  the  angle  of  attack.  The  pitching  moment  coefficient  lags  the 
angle  of  attack,  although  the  amount  of  lag  decreases  as  the  blade 
loading  increases.  These  trends  are  consistent  with  those  observed 
in  the  present  experimental  study  of  the  inlet  distortion  problem. 

The  similarities  exist  despite  the  differences  in  blade  profiles. 


j 


(Although  the  blades  in  the  cascade  test  have  a similar  thickness 
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distribution,  they  have  only  10*  of  camber  as  opposed  to  50*  for  the 
LSRR  blades.)  This  is  a strong  indication  that  cascading  effects  far 
outweigh  blade  profile  effects. 

Transient  Surface  Flow  Phenomena 

Hot-film  skin  friction  gages  were  mounted  adjacent  to  the 
pressure  transducers  to  detect  surface  flow  phenomena  such  as  separa- 
tion and  transition  to  turbulent  flow.  The  interpretation  of  the 
skin  friction  coefficient  derived  from  these  measurements  is  impos- 
sible without  some  reference  point.  This  reference  point  may  take 
the  form  of  a known  flow  condition  at  some  point  in  time  or  may  be 
based  on  the  results  from  a pressure  transducer  near  the  skin  friction 
gage  or  even  the  behavior  of  skin  friction  and/or  pressure  at  some 
other  point  on  the  blade.  For  example,  if  the  flow  is  known  to  be 
attached  and  laminar  at  some  point  in  time  and  the  skin  friction 
coefficient  subsequently  increases  and  shows  high  frequency  fluctua- 
tions, transition  to  turbulent  flow  is  indicated.  Conversely,  if  the 
skin  friction  drops  abruptly  and  shows  these  fluctuations,  separation 
of  the  boundary  layer  may  be  inferred.  If  the  flow  does  separate,  it 
may  reattach  downstream.  If  it  does  so,  the  flow  downstream  of  the 
reattachment  point  must  be  turbulent  and  will  be  characterized  by 
high  frequency  variations  in  the  skin  friction  coefficient  and  a 
generally  higher  average  value. 


i 


I 

I 

] 


i 


y 


86 

Most  of  the  activity  in  the  pressure  coefficient  and  skin 
friction  time  histories  is  confined  to  the  leading  edge  region; 
therefore,  in  this  discussion  the  pressure  coefficients  for  only  the 
first  two  measuring  stations  will  be  considered,  and  only  suction 
surface  data  are  presented.  Decreasing  pressure  coefficient  indicates 
increasing  normal  force.  Only  the  time-varying  part  of  the  signals 
are  plotted.  Figures  31  through  33  present  the  data  from  points  7. A, 
7.9,  and  7.15  respectively.  Figure  31  represents  the  lowest  average 
inlet  angle  tested.  This  results  in  the  highest  blade  loading  and 
results  in  separated  flow  during  part  of  the  cycle.  First  consider 
the  skin  friction  trace  for  x/c  * 0.03A.  At  the  point  marked  A,  the 
flow  is  attached  and  probably  laminar.  The  relatively  rapid  increase 
near  point  B may  be  taken  to  indicate  transition  to  turbulent  flow 
followed  almost  immediately  by  separation.  Separation  occurs  slightly 
before  minimum  inlet  angle.  (Recall  that  minimum  is  associated 
with  maximum  load.)  As  the  inlet  angle  increases  from  the  minimum 
value,  the  flow  reattaches  (point  C)  and  the  cycle  repeats  itself. 

It  is  interesting  to  note  that  the  separation  pattern  is  different 
for  each  cycle  plotted.  For  the  first  occurrence  of  separation  at 
the  extreme  left-hand  side  of  the  figure,  the  flow  at  x/c  ■ 0.034  is 
separated  over  a greater  range  of  inlet  angles  than  the  second 
occurrence.  The  third  time  the  flow  separates  (point  D)  it  does  not 
appear  to  reattach  after  minimum  inlet  angle  has  been  psssed.  The 
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pressure  transducer  at  this  chordwise  station  does  not  show  a signi- 
ficant decrease  in  suction  during  the  time  the  flow  is  separated,  but 
it  does  indicate  an  increase  in  the  magnitude  of  the  pressure  fluc- 
tuations. The  lack  of  a sharp  drop  in  leading  edge  suction  may  be 
taken  to  indicate  that  the  blade  has  made  only  a shallow  penetration 
of  the  dynamic  stall  regime.  That  is,  although  the  surface  flow  has 
become  separated,  the  external  potential  flow  is  relatively  unchanged. 
In  steady  flow,  the  blade  stalls  at  an  inlet  angle  of  28.5  degrees. 
Since  the  inlet  angle  in  distortion  drops  as  low  as  25*  it  is  evident 
that  the  blade  has  gone  past  the  quasi-steady  stall  point  by  a 
significant  amount.  Very  little  boundary  layer  activity  is  indicated 
by  the  skin  friction  gage  at  x/c  * 0.169.  During  the  time  that  the 
flow  at  x/c  * 0.034  is  separated,  there  is  slight  increase  in  the 
average  level  of  the  signal  and  in  the  amplitude  of  the  high  frequency 
fluctuations.  This  may  be  due  to  the  reattachment  of  the  flow 
between  x/c  * 0.034  and  x/c  * 0.169.  The  flow  downstream  of  a 
reattachment  point  is  expected  to  be  turbulent,  and  since  there  is  no 
indication  of  transition,  the  flow  must  always  be  turbulent  at  this 
stations.  The  pressure  transducer  trace  at  x/c  * 0.169  also  exhibits 
an  increase  in  the  high  frequency  fluctuations.  Taken  together,  all 
of  these  readings  indicate  that  a separation  bubble  probably  forms  on 
the  leading  edge  as  the  inlet  angle  approaches  a minimum  and  grows 
until  the  minimum  angle  has  been  passed  whereupon  it  shrinks  and  then 
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vanishes.  During  the  experiment  it  was  found  that  this  flow  condition 
results  in  the  lowest  average  inlet  angle  that  can  be  attained 
without  inducing  rotating  stall.  Thus,  in  the  case  of  inlet  distor- 
tion, the  onset  of  rotating  stall  appears  to  be  triggered  by  the 
occurrence  of  dynamic  stall  in  the  low  velocity  region  of  the  distor- 
tion. 

The  flow  condition  for  point  7.9  is  such  that  separation  does 
not  occur.  The  variations  in  the  skin  friction  trace  for  x/c  * 0.034 
shown  in  Fig.  32  are  due  to  transition  between  laminar  and  turbulent 
flow.  The  regions  of  higher  skin  friction  correspond  to  times  during 
which  the  transition  point  is  ahead  of  that  station  and  the  regions 
of  lower  skin  friction  correspond  to  times  during  which  the  transition 
point  is  behind  it.  When  the  flow  is  turbulent  at  x/c  = 0.034,  it 
must  also  be  turbulent  at  x/c  ■ 0.169.  Since  there  is  no  indication 
of  transition  at  x/c  = 0.169,  the  flow  must  always  be  turbulent 
there.  Thus,  the  transition  point  must  always  lie  ahead  of  x/c  = 

0.169.  At  the  lower  values  of  instantaneous  inlet  angle,  the  transi- 
tion point  is  ahead  of  x/c  * 0.034.  The  data  obtained  from  point 
7.15  are  shown  in  Fig.  33.  The  small  humps  in  the  skin  friction 
trace  for  x/c  ■ 0.034  indicate  transition  to  turbulent  flow  only  at 
the  minimum  inlet  angle.  By  the  same  reasoning  as  before,  the  flow 


at  x/c  ■ 0.169  must  always  be  turbulent. 
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TABLE  IV 

HARMONIC  CONTENT  OF  UNSTEADY  ROTOR  PRESSURE  RISE  COEFFICIENT 


a0 

al 

♦l 

a2 

♦ 2 

Point 
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.37.1 

0.013 

70.9 
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0.090 
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0.017 

335.6 
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0.393 
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0.010 

17.0 

a3 

♦3 
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pressure,  normal  force,  and  pitching  moment  coefficient  loops.  There 
is  moderate  second  harmonic  content  in  the  loops  for  points  7.4  and 


7.9.  However,  the  loop  for  point  7.15,  which  is  in  the  potential 
flow  regime  and  also  in  the  low  loss  region,  is  almost  a pure  first 
harmonic . 

Velocities  and  Flow  Angles  at  the  Rotor  Exit  Plane 

Three  velocity  components  were  measured  approximately  0.07 
chord  lengths  aft  of  the  rotor  exit  plane  using  a triaxial  hot-film 
probe.  The  parameters  which  will  be  employed  to  present  the  results 
of  these  measurements  will  be  the  exit  velocity,  Wj , the  circumfer- 
ential flow  angle,  82,  and  the  radial  flow  angle,  We  will  first 

consider  measurements  made  in  undistorted  flow  at  an  inlet  angle  of 
42.3*.  The  variation  in  the  exit  velocity  across  the  blade  gap  is 
plotted  in  Fig.  35.  The  point  y/g  ■ 0.0  corresponds  to  the  downstream 
projection  of  the  suction  surface  of  one  blade  and  y/g  “ 1.0  corre- 
sponds to  the  downstream  projection  of  the  pressure  surface  of  its 
neighbor.  The  velocity  profile  is  relatively  flat  showing  a slight 
decrease  as  one  approaches  the  pressure  surface  (y/g  * 1.0).  There 
is  an  abrupt  drop  in  velocity  near  the  blade  trailing  edges  at  y/g  “ 
0.0  and  y/g  ■ 1.0.  It  must  be  pointed  out  that  the  sensor  arrangement 
of  this  particular  probe  is  such  that  is  is  not  suitable  for  use  in 
regions  where  strong  velocity  gradients  exist.  The  blade  wakes  are 


obviously  such  regions;  however,  the  exit  velocity  should  be  only 


97 

slightly  affected.  The  radial  flow  angle  will  be  affected  somewhat 
more  and  the  circumferential  flow  angle  will  be  strongly  affected. 

The  circumferential  exit  angle  (Fig.  36)  is  constant  except  in  the 
blade  wake  regions.  The  flow  angle  in  these  regions  are  probably 
erroneous  due  to  the  previously  mentioned  limitations  of  the  probe. 

The  radial  flow  angle  (Fig.  37)  is  extremely  close  to  zero  except  in 
the  blade  wake  regions. 

In  order  to  study  the  velocity  and  flow  angles  in  the  case  of 
distorted  flow,  we  will  again  consider  the  four  inlet  angles  corres- 
ponding to  maximum  axial  velocity,  minimum  axial  velocity  and  the 
average  axial  velocity  with  the  inlet  angle  increasing  and  with  it 
decreasing.  The  velocity  profile  for  these  angles  is  plotted  in  Fig. 

38.  Little  variation  in  the  exit  velocity  is  noted  because  it  is 
composed  of  both  an  axial  part  (which  is  varying)  and  a circumferen- 
tial part  due  to  the  rotation  of  the  rotor  (which  is  relatively 
constant).  One  notes  that  the  two  cases  which  have  the  same  inlet 
angle  do  not  have  identical  velocity  profiles.  The  plots  indicate 
that  the  velocity  distortion  has  been  displaced  in  the  direction  of 
rotation  of  the  rotor.  The  variation  in  Bj  across  the  gap  is  shown 
in  Fig.  39.  This  profile  is  relatively  flat  except  for  an  inlet 
angle  of  43.8*  with  the  inlet  angle  decreasing.  In  this  case,  the 
exit  angle  exhibits  rather  pronounced  fluctuations  for  which  there  is 

I 

presently  no  explanation.  It  is  apparent  that  the  gap-averaged 
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circumferential  exit  angle  is  undergoing  significant  changes  as  the 
distortion  is  traversed.  By  obtaining  gap-averaged  values  of  the 
exit  angle  at  these  and  several  additional  locations  in  the  distortion 
field  and  plotting  them  as  a function  of  the  instantaneous  inlet 
angle  (Fig.  40)  one  generates  a loop  similar  to  those  produced  from 
the  pressure  data.  Comparing  this  loop  to  the  undistorted  exit  angle 
behavior,  one  concludes  that  the  unsteady  exit  angle  is  lagging  the 
steady  curve.  The  shape  of  the  loop,  however,  is  a very  poor  match 
to  the  shape  of  the  steady  state  curve.  Furthermore,  the  variations 
in  exit  angle  measured  by  the  probe  appear  to  be  excessive.  Exit 
angles  were  also  calculated  by  a momentum  theory  approach  (similar  to 
that  of  Appendix  V)  neglecting  any  unsteady  effects  not  inherent  in 
the  normal  force  and  rotor  pressure  rise  loops.  The  results  of  these 
calculations  (which  are  not  presented  here)  tend  to  support  the 
belief  that  the  indicated  exit  angle  variations  in  distorted  flow  are 
excessive.  The  radial  flow  angle  (Fig.  41)  is  definitely  not  zero 
indicating  that  the  distorted  flow  field  has  induced  three-dimensional 
effects  into  the  rotor  wake.  The  lack  of  agreement  between  the  exit 
angle,  @2>  measured  in  distorted  flow  and  that  in  uniform  inlet 
flow  may  be  partly  due  to  these  three-dimensional  effects. 


/3,  - 46.9° 
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CONCLUSIONS 

An  inlet  distortion  experiment  has  been  performed  in  which 
an  isolated  compressor  rotor  was  subjected  to  a sinusoidal  variation 
in  inlet  axial  velocity.  The  instrumentation  system  developed  for 
this  experiment  allows  the  measurement  of  unsteady  blade  pressures  on 
the  rotor  with  a high  degree  of  precision  and  without  the  noise 
problems  normally  associated  with  on-rotor  measurements.  Qualitative 
skin  friction  measurements  were  also  performed  on  the  rotor  blading. 

A wake  probe  mounted  in  the  rotating  reference  frame  was  used  to  make 
velocity  and  flow  angle  measurements  at  the  rotor  exit  plane.  The 
development  of  the  capability  to  make  such  on-rotor  measurements  with 
a precision  and  noise-immunity  commonly  possible  only  for  non-rotating 
experiments  is  a significant  accomplishment  of  the  present  program. 
From  the  experimental  results  which  have  been  presented  herein,  the 
following  conclusions  may  be  drawn. 

1.  The  peak  pressure  rise  on  a circumferentially  averaged  basis 
is  greater  for  undistorted  flow  than  for  distorted  flow. 
However,  at  low  average  inlet  angle  (high  blade  loading), 
the  rotor  produces  a slightly  higher  pressure  rise  in 
distorted  flow.  This  is  attributed  to  the  fact  that  the 
blade  transiently  operates  beyond  the  steady  state  stall 
point  at  these  average  inlet  angles  and  yet  does  not  stall. 
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2.  Under  high  blade  loading, 'the  blade  pressure  distribution 
exhibits  a sharp  negative  'peak  on  the  suction  surface  near 
the  leading  edge.  This  peaking  is  more  pronounced  than  is 
normally  the  case  with  an  isolated  airfoil. 

3.  As  the  inlet  angle  is  increased,  resulting  in  reduced  blade 
loading,  the  pressures  on  the  rear  portion  of  the  blade 
remain  essentially  unchanged.  On  the  leading  edge,  however, 
the  pressures  undergo  considerable  change.  At  the  very  low 
blade  loadings  the  pressures  on  the  suction  surface  of  the 
leading  edge  are  greater  than  those  on  the  pressure  surface 
resulting  in  negative  loading  on  the  leading  edge. 

4.  As  the  blade  traverses  the  distortion,  the  greatest  pressure 
variations  occur  at  the  most  forward  measuring  station  (x/c 
* 0.034).  The  response  at  this  station  consistently  lags 


the  inlet  angle.  The  response  at  the  second  measuring 
station  (x/c  * 0.169)  is  much  reduced  in  amplitude  and  leads 


the  inlet  angle  at  high  blade  loading,  changing  to  a lagging 
response  as  the  blade  loading  decreases.  Pressures  at  the 
third  and  fourth  measuring  stations  (x/c  • 0.381  and  x/c  ■ 
0.619,  respectively)  show  very  little  variation  for  quasi- 
steady flow.  However,  under  the  influence  of  the  unsteadi- 
ness created  by  the  inlet  distortion,  a moderate  response  is 
observed.  The  pressures  farther  back  on  the  blade  essentially 
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follow  the  quasi-steady  characteristic  and  exhibit  very 
little  change  with  varying  inlet  angle. 

5.  The  normal  force  coefficient  loops  show  a significant 
departure  from  the  quasi-steady  characteristic.  The  ampli- 
tude of  the  response  is  increased  considerably  and  it 
generally  leads  the  inlet  angle.  A significant  penetration 
beyond  the  steady-state  stall  point  is  indicated. 

6.  The  pitching  moment  coefficient  loops  are  dominated  by  the 
pressure  response  at  the  leading  edge.  They  show  good 

agreement  with  the  steady  state  curve  at  high  inlet  angles 

\ 

but  depart  somewhat  from  the  steady  state  at  low  inlet 
angles . 

7.  For  the  lowest  average  inlet  angle  which  could  be  attained 
without  the  occurrence  of  rotating  stall,  a separation 
bubble  is  observed  to  exist  on  the  leading  edge  of  the  blade 
while  it  is  in  the  low  velocity  (high  loading)  region.  This 
may  indicate  that  the  inception  of  rotating  stall  is  trig- 
gered by  the  occurrence  dynamic  stall  in  the  low  velocity 
region. 

8.  Velocity  and  flow  angle  measurements  made  at  the  rotor 
exit  plane  at  blade  midspan  indicate  that  the  flow  is 
essentially  two-dimensional  in  the  undistorted  case.  The 
wake  profiles  are  smooth  and  relatively  constant  except  for 


I 
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the  regions  immediately  behind  the  blades  where  there  is  a 
velocity  deficit  and  some  flow  angle  variations.  In  distor- 
ted flow,  however,  significant  radial  flow  is  indicated. 

This  radial  velocity  component  changes  direction  (from 
inward  to  outward)  as  the  distortion  is  traversed.  It  thus 
appears  unreasonable  to  attempt  to  make  quantitative  measure- 
ments downstream  of  a rotor  undergoing  distortion  without 
performing  a radial  survey  of  the  flow  parameters. 

9.  The  circumferential  exit  angle  variations  measured  in 

distorted  flow  appear  to  be  excessive.  It  is  possible  that 
an  interaction  between  the  probe  and  the  unsteady  flow  field 
behind  the  rotor  is  occurring. 

The  following  additional  statements  may  be  made  concerning  the 
instrumentation  system. 

10.  The  use  of  on-board  signal-conditioning  and  amplifying 
equipment  produces  signals  which  are  not  contaminated  by 
slip-ring  noise  and  reduces  the  number  of  slip-rings  required 
for  each  channel  of  data. 

11.  The  circuitry  which  was  developed  to  perform  these  signal- 
conditioning and  amplifying  functions  proved  to  be  extremely 
reliable. 

12.  Tnc  triaxial  hot-film  probe  usqd  to  measure  flow  angels  in 
distorted  flow  was  not  satisfactory.  The  probe  geome  try  (in 
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particular,  the  sensor  placement)  is  believed  to  be  at 
fault.  In  any  case,  a better  technique  for  measuring  these 
flow  angles  is  needed. 


no 


APPENDIX  I 
POWER  SUPPLIES 

The  unregulated  D.  C.  supply  for  the  rotating  frame 
instrumentation  was  located  outside  the  LSRR.  Two  20  volt  supplies 
capable  of  delivering  a current  of  3 amperes  were  connected  as 
indicated  in  Fig.  42  to  provide  plus  and  minus  twenty  volts  with 


FIG.  42  UNREGULATED  D.C.  POWER  SUPPLY 

respect  to  a common  ground.  Two  types  of  voltage  regulators  were 
used  in  the  rotating  frame  ins trument at  ion  package.  The  first  type 
was  a dual  regulator  supplying  + 13  volts  to  the  pressure  transducer 
amplifiers  and  the  constant  temperature  anemometers.  Its  schematic 
diagram  is  shown  in  Fig.  43.  This  circuit  was  taken  directly  from 


Ref.  37. 
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FIG.  43  SCHEMATIC  DIAGRAM  OF  ± 15  VOLT  REGULATOR 

The  output  capacitors,  C2  and  C4  are  required  for  stability  and 
to  reduce  high  frequency  noise.  The  input  capacitors,  Cl  and  C3 , are 
necessary  because  the  regulator  is  located  some  distance  from  the 
filter  capacitors  in  the  unregulated  supply.  Capacitors  Cl  and  C2 
are  ceramic  disk  types  while  C3  and  C4  must  be  solid  tantalum.  The 
diodes  D1  and  D2  allow  the  regulators  to  start  up  under  a common  load 
(i.e.,  one  between  VQ*t  and  VQ~t)  when  the  unregulated  supply 
is  turned  on. 

The  circuit  for  the  pressure  transducer  excitation  supply  of  -10 
volts  was  also  taken  from  Ref.  37.  Its  schematic  diagram  appears  in 


Fig.  44. 
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| Capacitors  Cl  and  C2  are  solid  tantalum  with  Cl  being  required 

because  the  regulator  is  remote  from  the  unregulated  supply  and  C2 
being  required  for  stability.  The  resistors  R1  and  R2  form  a divider 
which  causes  the  output  to  be  twice  the  basic  regulator  output  of  -5 
volts.  The  bypass  capacitor  C3  improves  ripple  rejection  and  transient 
^ response.  It  may  be  an  aluminum  electrolytic. 
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APPENDIX  II 

PRESSURE  TRANSDUCER  AMPLIFIERS 
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In  order  to  boost  the  output  of  the  pressure  transducers  mounted 
on  a rotor  blade  to  a level  where  it  would  not  be  susceptible  to 
noise  contamination,  it  was  necessary  to  place  a fixed  gain  amplifier 
inside  the  rotor  hub  for  each  transducer.  A gain  of  100  was  chosen 
based  on  the  expected  pressure  levels  and  transducer  sensitivities  as 
well  as  the  restrictions  which  the  stationary  frame  equipment  placed 
on  the  output  voltage  (no  more  than  +_  1 .A  volt  to  avoid  overloads). 

The  circuit  which  was  developed  to  serve  this  function  is  shown  in 
Fig.  45.  Capacitors  C4  through  C7  are  merely  filters  used  to  minimize 
ripple  on  the  operational  amplifiers  supply  lines.  Amplifiers  A1  and 
A2  are  connected  as  unity  gain  voltage  followers  whose  sole  purpose 
is  to  present  a high  input  impedance  to  the  transducer  bridge. 
Amplifier  A3  and  resistors  R1  through  R6  form  a differential  amplifier 
with  a gain  of  100.  The  fixed  resistors  are  1 percent  metal  film 
types  and  the  variable  resistors  are  22-turn  cermet  trimmers.  These 
types  were  chosen  for  their  stability  with  age  and  temperature.  The 
output  of  the  amplifier  may  be  shifted  by  adjusting  R6.  Since  each 
amplifier  is  always  connected  to  the  same  transducer,  this  adjustment 
is  done  only  once  during  the  initial  setup  to  account  for  any  offset 


in  the  transducer.  Trimmer  R4  is  adjusted  to  peak  the  D.  C.  common 
mode  rejection.  The  offset  adjustment  (R6)  and  the  D.  C.  common  mode 
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adjustment  (R4)  interact  slightly  but  this  was  not  considered  a 
problem  because  they  are  set  only  once.  The  D.  C.  common  mode 
rejection  is  considered  to  be  properly  set  when  a D.  C.  voltage 
applied  simultaneously  to  inputs  A and  B does  not  change  the  output. 
After  these  adjustments  are  completed,  the  AC  common  mode  is  peaked. 
Capacitors  Cl  and  C2  are  nominally  20  pf.  Additional  capacitors  (on 
the  order  of  a few  picofarads)  are  paralleled  with  one  or  the  other 
to  maximize  the  AC  common  mode  rejection.  A function  generator  was 
used  to  apply  a sine  wave  simultaneously  to  both  inputs  and  a trial 
and  error  process  was  followed  until  the  output  ripple  was  minimized 


where  G is  nominally  100.  Since  each  amplifier  is  always  connected 
to  the  same  transducer,  the  exact  value  of  G is  unimportant.  Rather, 
an  output  sensitivity  of  the  transducer-amplifier  pair  is  experi- 
mentally determined  and  this  sensitivity  is  used  in  all  calculations. 

Reference  38  is  an  exhaustive  treatment  of  the  theory  and 
applications  of  operational  amplifiers  which  was  written  for  the 
layman.  It  is  recommended  reading  for  anyone  contemplating  a similar 
use  of  operational  amplifiers. 
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APPENDIX  III 

OONSTANT  TEMPERATURE  ANEMOMETERS 

The  constant  temperature  anemometer  units  used  in  this  experiment 
are  an  original  design  initially  developed  for  use  with  the  skin 
friction  gages.  After  the  commercial  units  which  had  been  purchased 
to  operate  the  triaxial  hot-film  probe  failed  under  centrifugal 
acceleration  loading,  the  in-house  design  was  used  for  that  applica- 
tion also.  The  units  were  designed  to  be  inexpensive,  rugged,  and 
stable  because  their  location  inside  the  hub  would  subject  them  to 
centrifugal  loads  and  would  preclude  easy  adjustment  once  installed. 
These  goals  were  easily  met  by  the  circuit  of  Fig.  46.  First  consider 
the  bridge  composed  of  R12  through  R15.  When  the  voltage  at  point  A 
is  the  same  as  the  voltage  at  point  B the  bridge  is  said  to  be 
balanced.  The  output  of  the  error  amplifier  composed  of  A1 , R6 , R7 , 
R16,  and  R17  will  be  aero.  The  output  of  A2  will  be  controlled  by 
the  setting  of  R1  which  allows  the  circuit  to  have  an  output  (at 
point  C)  with  a aero  error  signal.  If  this  feature  were  emitted,  the 
anemometer  would  not  start  up  properly  when  power  was  applied.  The 
transistors  Q1  and  Q2  connected  in  a Darlington  configuration  boost 
the  current  output  of  A2  to  meet  the  demands  of  the  bridge.  The 
Darlington  pair  is  stabilised  by  resistor  R1 1 . Returning  to  the 
bridge,  one  finds  that  the  condition  for  balance  is  that 


FIG.  46  SCHEMATIC  DIAGRAM  OF  CONSTANT  TEMPERATURE  ANEMOMETER 
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R12  _ R13. 
R14  = Rl 5 


(35) 


or  alternately 


R15  _ RJL3. 
R14  = R12 


(36) 


Thus  it  can  be  seen  that  the  bridge  is  balanced  when  the  probe  resis- 
tance ( Rl 5 ) is  one  fifth  that  of  the  control  resistor  (R14).  (This 
ratio  is  set  by  the  chosen  values  of  R13  and  R12  and  is  what  is 
commonly  referred  to  as  a 5:1  bridge.)  The  probe  is  being  heated  by 
the  current  through  it  and  is  operating  at  a resistance  somewhat 
higher  than  its  cold  resistance.  If  for  some  reason  it  were  to 
become  cooler  than  the  temperature  at  which  the  bridge  balanced,  its 
resistance  would  drop  causing  the  voltage  at  point  A to  be  lower  than 
that  at  poiint  B.  The  error  amplifier  would  then  develop  a negative 
output  which  would  be  inverted  and  further  amplified  by  A2  resulting 
in  a higher  voltage  at  point  C.  The  increased  current  through  the 
probe  would  heat  it  causing  it  to  return  to  a balanced  condition. 
Actually,  the  bridge  will  not  always  be  perfectly  balance;  however, 
the  imbalance  will  be  inversely  proportional  to  the  system  gain. 

This  gain  may  be  adjusted  over  a range  of  roughly  50  to  1000  by  R8  so 
for  all  practical  purposes  the  sensor  will  operate  at  the  desired 
temperature.  A more  detailed  description  of  the  operation  of  constant 
temperature  anemometers  may  be  found  in  Ref.  39.  A brief  discussion 
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of  some  practical  aspects  of  anemometry  and  an  extensive  bibliography 
of  related  articles  may  be  found  in  Ref.  40. 
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APPENDIX  IV 

INSTRUMENTATION  USED  AND  FLOW  CONDITIONS  RUN 

Table  V details  the  instrumentation  available  in  configurations 
one  and  two.  Note  that  the  stationary  frame  instrumentation  is  the 
same  for  both  configurations.  Configuration  one  was  used  for  runs 
six  through  eight  and  configuration  two,  for  runs  ten  and  eleven. 

The  hot  film  probes  behind  the  rotor  were  used  to  measure  the  number 
and  velocity  of  stall  cells  when  rotating  stall  was  occuring.  The 
thermocouple  was  intended  to  detect  any  rapid  temperature  variations 
which  would  not  be  sensed  by  the  static  data  system  so  that  the  point 
could  be  repeated.  Table  VI  is  a listing  of  the  flow  conditions  run 
in  the  present  experimental  program.  The  flow  condition  is  labelled 
as  undistorted  (U)  or  distorted  (D).  Judgement  as  to  the  existence 
of  rotating  stall  was  based  on  observations  of  the  unsteady  instru- 


mentation output  on  an  oscilloscope. 
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TABLE  V 

INSTRUMENTATION  CONFIGURATION  FOR  LSRR 


Rotating  Frame 

Instrument  at  ion 

Channel 

Conf.  1 

Conf.  2 

1 

XDl 

XDl 

2 

XD2 

XD2 

3 

XD3 

XD3 

4 

XD4 

XD4 

5 

XD5 

XD5 

6 

XD6 

XD6 

7 

XD7 

SF1 

8 

XD8 

SF2 

9 

XD9 

SF3 

10 

XD10 

SF4 

11 

XDl  1 

SF5 

12 

SF1 

SF6 

13 

SF2 

TX1 

14 

SF3 

TX2 

13 

SF5 

TX3 

16 

SF8 

PT 

17 

SF1 1 

XDl  2 

Stationary  Frame 

Instrument  at  ion 

Channel 

Data 

1 

HF1 

2 

HF2 

3 

HF3 

4 

OPR 

5 

TC 
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TABLE  V (Cont'd) 

INSTRUMENTATION  CONFIGURATION  FOR  LSRR 


Key:  XD 

SF 
TX 
PT 
HF 
OPR 
TC 


Pressure  Transducer 

Skin  Friction  Gage 

Triax  Probe  Sensor 

Total  Pressure  in  Wake 

Hot  Film  Probe  Behind  Rotor 

Once /Revolution 

Thermocouple  Behind  Rotor 
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TABLE  VI 

SUMMARY  OF  FLOW  CONDITION  USED  IN 
LSRR  INLET  DISTORTION  EXPERIMENT 


Kun 

Pc'  i nt 

• i1 

® MCI. 

1* 

W 

l«/l. 

Flow 

Condi t ion 

( degf  •»  r*  sc  ) 

< degrees  ) 

( f t /sec  ) 

( f t / SPC  ) 

( f t / sec  ) 

6 

b 

2b  . ** 

1.9 

78 .9 

144 . 7 

164.6 

0.34  3 

U 

6 

7 

33.2 

-3.2 

94  9 

143. 1 

173.4 

0.6  34 

t’ 

6 

in 

3b .«» 

-6.0 

103.3 

144.9 

179.1 

0.727 

l’ 

6 

1 1 

38  . 6 

-8.8 

116.2 

144 . 7 

183.6 

0.803 

V 

6 

i: 

4() . 3 

-10.3 

122.8 

149.0 

190.0 

0.847 

V 

6 

13 

^3.  1 

-13.1 

139.9 

143.  1 

198.8. 

0.937 

l 

b 

18 

4m  . 7 

-1m. 7 

143. 1 

144.8 

203.6 

0.988 

l' 

b 

2) 

m7 .9 

-17.9 

160.3 

143.0 

216.3 

J .107 

u 

b 

24 

V>.2 

-2U.2 

173.8 

144  .8 

226.2 

1.200 

u 

7 

28.9 

1 . 1 

78.1 

141.7 

lb  1 .8 

0.331 

*D 

7 

33.9 

-3.3 

94  .0 

1 4 1 .8 

170.1 

0.663 

D 

7 

9 

3t>.4 

-b  . 9 

106  .b 

M .8 

177.4 

0.732 

D 

7 

i: 

40  . m 

-10.4 

120.4 

141  .7 

183.9 

0.8  30 

D 

7 

42 . 2 

-12.2 

128.3 

141 .7 

191  .1 

0.903 

D 

7 

16 

43.7 

-13.7 

* 133.9 

U1  .8 

196.1 

0.936 

D 

7 

19 

49.2 

-13.2 

142. 3 

141  .7 

201.0 

1 .006 

D 

7 

22 

28.7 

1.3 

77.9 

141  .7 

161  .3 

0.347 

o, 

St  al  led 

7 

23 

2b  .u 

3.6 

70.3 

141 .7 

138.2 

0.49b 

o, 

St  si  led 

7 

24 

30.7 

-0.7 

(U.  . 1 

141.7 

164.8 

0.394 

D, 

Seal  led 

a 

3 

28 .8 

1.2 

79 . 3 

144  .0 

164  .4 

0.331 

u, 

St  at ltd2 

H 

4 

29.8 

0.? 

82.3 

143.9 

163.8 

0.372 

u, 

Stal ltd2 

H 

s 

30.  b 

-o . a 

83.0 

143.9 

167.1 

0.391 

l\ 

Seal  ltd2 

in 

3-14 

43.  * 

-13.8 

1 36  .a 

142.6 

197.3 

0.938 

D3 

Stal  ltd3 

10 

18-2* 

27.9 

2.1 

73.3 

142.4 

161 . 1 

0.329 

D, 

1 1 

3- la 

42.3 

-12.3 

129.0 

142.0 

191  .8 

0.908 

u3 

1 1 

19-31 

31  .3 

-1.9 

87.1 

141  .9 

166.3 

0.614 

u* 

Stal ltd3 

1 1 

34-38 

32.0 

-2.0 

88.  H 

141  .9 

167.4 

0.626 

Not  e* 

( 1 ) 

(21 

6 1 is  based  on  ?R  for  distorted  flow 

After  coai^letion  of  run  it  was  discovered 

t hat  a 

bear  i ng 

had  overheated  and 

failed  melting  lose  of  the  plastic  tubing 

running 

to  the 

(3) 

Wake  probe  traverses  at  a single  flow  condition 
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APPENDIX  V 

CALCULATION  OF  EXIT  ANGLE  USING  MOMENTUM  THEORY 

One  may  derive  an  expression  relating  the  exit  angle  to  the 
normal  force  coefficient  of  the  blades,  the  static  pressure  rise 
across  the  rotor,  and  the  inlet  angle  through  momentum  considerations. 

The  control  volume  which  will  be  used  in  the  subsequent  deriva- 
tion is  shown  in  Fig.  47.  Note  that  the  streamline  contraction  due 
to  boundary  layer  growth  is  exaggerated.  The  continuity  equation  may 
be  written  in  terms  of  axial  velocities, 

PlcXlSbl  = p2Cx2gb2  (37) 

and  for  incompressible  flow,  pj  “ p2  “ p,  and 

CXlbl  ” CX2b2  • (38) 

The  tangential  momentum  equation  is 

Fe  “ pg  (cXlce1bi  ~ Cx2ce2b2)  (39) 

which  may  be  simplified  using  Eq.  (38)  to  read 

F0  " PRCx,bl(C0j  ~ C02)  (40) 

The  axial  momentum  equation  requires  some  assumption  as  to  the 
variation  in  pressure  through  the  rotor.  If  one  assumes  a linear 
variation,  the  axial  momentum  equation  becomes 

P8  + P8 

Fx  * Ps,8b2  - Ps,8bl  + — 8<bl  ~ b2> 

2 

+ pgb2Cy2  ~ PgbiCyj  . 


(41) 
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The  first  term  on  the  right  side  of  the  equality  is  due  to  pressure 
on  the  rear  face  of  the  control  volume,  the  second  is  due  to  pressure 
on  the  front  face,  and  the  third  is  due  to  pressure  on  the  sides. 

The  last  two  terms  represent  the  axial  momentum  change  due  to  the 
velocity  difference.  Simplifying, 


g( 


b 1+  b2 


■)(P8. 


- PSl)  + P gb l C: 


X1 (Cx2 


- cXl) 


(42) 


The  axial  and  tangential  forces  may  be  combined  to  yield  the  normal 
force  by  using 

Fn  = Fe  sin£  + Fx  cose  (43) 

where  e is  the  chordal  stagger  angle.  The  normal  force  coefficient 
is  defined  as 

« -Fjj_  (44) 


Inserting  Eqs.  (40)  and  (42)  into  Eq.  (43)  and  using  the  definition 
of  Cjg  as  given  by  Eq.  (44)  one  obtains 


CN 


g 

4(-)sin2B  sine  cotB,  - cotB  (- — ) 

cl  1 *■  ^2 


This  may  be  solved  for  82  as 


127 


cote. 


/ x2 

C CN^  + r1) 

*2  „ Lx, 

cT°t6i  ~ ~ 

^1  z./^ 


4(^)sin261sinC 


(46) 


, x2 
V1  * ^ 

4 tan26i  tan^  'C^  ^ 


+ (Jli)(J!2.  1)cotC 


This  equation  was  used  to  obtain  the  exit  flow  angle,  82 , from 
measured  values  of  inlet  angle,  normal  force  coefficient,  and  rotor 
pressure  rise  coefficient. 


I 
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APPENDIX  VI 

TRIAXIAL  PROBE  CALIBRATION 

The  probe  used  in  the  LSRR  inlet  distortion  experiment  is  a 
Thermo-Systems,  Inc.  1296M  triaxial  hot-film  probe.  A total  pressure 
transducer  is  attached  to  the  support  shaft.  Two  axis  systems  are 
used  in  calibrating  the  triaxial  hot-film  probe.  One  system  is  the 
axial-radial-circumferential  axis  system  of  the  rotor  while  the  other 
is  chosen  with  the  three  axis  parallel  to  the  three  mutually  ortho- 
gonal velocity  sensors  which  make  up  the  probe.  These  two  systems 
are  shown  in  Fig.  48  omitting  the  total  pressure  probe  for  clarity. 


vi 
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FIG.  48  AXIS  SYSTEMS  USED  IN  WAKE  PROBE  CALIBRATIONS 


12) 

The  velocity  components  in  the  two  systems  may  be  related  by  two 
systems  may  be  related  by  two  sets  of  transformation  equations.  The 
first  set  gives  velocities  in  the  rotor  reference  frame,  v^,  in 


terms  of  the  sensor  frame  velocities,  u^ . 

Vi  = anU) 

+ a2 iu2  + a3 1U3 

(47) 

v2  « ai2ui 

♦ a22u2  + a32U3 

(48) 

v3  ♦ anui 

+ a23u2  + a33u3 

(49) 

The 

second  set  of  equations  is  the  reverse  transformation 

Ui  = a i i Vi 

♦ ai 2v2  + ai  3V3 

9 

(50) 

u2  + a2 i Vi 

♦ a22v2  + a2  jv 3 

i 

(51) 

and 

u 3 = anVi 

+ a 3 2 v2  ♦ a33V3 

(52) 

The 

coefficients,  a^j, 

used  in  these  equations  are  as 

fol lows . 

aii  = /277 

a 2 1 = - /1/6 

a3i  = 

- /I/ 6 

ai2  = 0 

a2 2 * - /1/2 

a 3 2 = 

/T7T 

a 1 3 = /!/ 3 

a2  3 = /!/$ 

a 3 3 = 

/T7T 

The  probe  was  calibrated  in  the  NCSU  Low  Speed  Wind  Tunnel  using 


a support  fixture  which  allowed  the  probe  to  be  pitched  and  yawed 
while  maintaining  the  sensors  in  precisely  the  same  location  within 
the  tunnel,  thus  allowing  any  flow  nonuniformities  which  might  exist 
in  the  tunnel  to  be  neglected.  The  sign  conventions  used  for  the 
pitch  and  yaw  angles,  e and  $,  respectively,  are  shown  in  Fig.  49. 
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FIG.  49  SIGN  CONVENTIONS  FOR  PITCH  AND  YAW  ANGLES 


These  angles  and  the  total  velocity  vector,  V,  may  be  used  to  calculate 
the  velocities  in  the  rotor  reference  frame. 


Vj  - V sin0 
v2  - V cos0  sin4> 


(53) 

(54) 


v 


V COS0  COSlfi 


(55) 
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The  direction  cosines  of  the  sensors  were  used  to  calculate 
pitch  and  yaw  angles  which  would  place  each  sensor  in  turn  perpendi- 
cular to  the  flow.  Based  on  previous  experience,  it  was  suspected 
that  the  probe  shaft  was  deflecting  the  flow  so  that  the  flow  was  not 
actually  normal  to  the  sensor  in  question  at  the  calculated  pitch  and 
yaw  angle.  Consequently,  with  the  tunnel  running,  the  probe  was 
yawed  slightly  while  monitoring  the  anemometer  output  voltage.  The 
flow  was  taken  to  be  normal  to  the  sensor  at  maximum  output  voltage. 
The  yaw  angle  at  peak  output  did  vary  slightly  from  the  calculated 
values  confirming  that  the  probe  shaft  was  indeed  interfering  with 
the  airflow  over  the  sensors.  With  each  sensor  in  turn  placed  normal 
to  the  flow,  the  tunnel  was  operated  over  a wide  range  of  velocities 
and  the  various  flow  properties  and  the  output  voltage  of  the  appro- 
priate anemometer  were  recorded.  It  was  desired  to  fit  these  data  to 
the  form  of  Eq.  (56). 

F.2  - El  « K/pV  (Ts  - TJ  , (56) 

where  E is  the  anemometer  output  voltage,  p is  the  density,  and 
Tg  and  T„  are  the  sensor  and  free  stream  temperature,  respectively. 
The  constants  EQ  and  K are  determined  from  the  calibration  data. 

This  relation  between  the  fluid  properties  and  the  anemometer  voltage 
was  first  suggested  by  King  (Ref.  41)  and  bears  the  name  King's  Law. 
This  type  of  relation  and  some  of  its  limitations  are  discussed  in 
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Refs.  40  and  42.  Once  such  limitation  is  that  no  single  set  of 
values  for  EQ  and  K is  valid  over  the  entire  velocity  range; 
consequently,  a least  squares  procedure  was  used  to  obtain  two  curve 
fits  of  the  form  of  Eq.  (56),  one  of  which  was  for  the  higher  veloci- 
ties and  the  other,  for  the  lower  velocities.  The  anemometer  voltage 
at  the  intersection  of  these  two  curves  was  calculated.  During  the 
remainder  of  the  calibration  procedure,  if  the  output  voltage  of  an 
anemometer  was  higher  than  the  crossover  voltage  the  high  velocity 
fit  for  that  sensor  was  used,  and  if  it  was  lower,  the  low  velocity 
fit  was  used.  Typical  data  obtained  for  one  of  the  sensors  and  the 
two  curve  fits  used  are  presented  in  Fig.  50. 

After  determining  the  response  of  the  probe  sensors  to  normal 
flow,  it  was  necessary  to  investigate  the  behavior  of  the  sensors  in 
non-perpendicular  flows.  One  expression  which  is  commonly  used  to 
describe  the  response  of  a sensor  to  transverse  flow  is: 


^ef  f 


u2  + K2u2 
n t 


(57) 


where  qeff  is  the  effective  velocity  measured  by  the  sensor,  un 
and  ut  are  the  normal  and  transverse  velocity  components,  respec- 
tively, and  K is  a constant  which  is  empirically  determined.  Unfor- 
tunately, this  expression  does  not  lend  itself  to  situations  where 
the  flow  may  be  blocked  by  the  probe  shaft  and  the  behavior  of  the 

4 

probe  indicated  that  there  were  serious  blockage  effects.  (See  Ref.  43 
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FIG.  50  TYPICAL  RESULTS  OF  NORMAL  VELOCITY  CALIBRATION 


78-1 1-140-6 
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for  a discussion  of  probe  shaft  interference  effects.)  The  formulation 
of  Eq.  (58)  was  developed  to  include  these  blockage  effects  as  well 
as  sensitivity  to  transverse  flow. 

q|ff  “ » (58) 

where  qeff  and  up  are  as  before  and  f(9,  ♦)  is  a weighting  factor 
of  order  unity  which  is  a function  of  the  pitch  and  yaw  angles.  The 
wind  tunnel  was  operated  at  a constant  velocity  of  116  ft/sec  while 
the  probe  was  pitched  and  yawed  and  the  data  obtained  used  to  calcu- 
late the  value  of  this  function  for  each  sensor  at  each  combination 
of  pitch  and  yaw  angles. 

The  values  of  the  function  at  each  value  of  9 were  fitted  using 
a least-squares  method,  to  a fifth  order  polynomial  in  4>.  The 
coefficients  of  each  term  in  this  polynomial  were  then  fitted  by 
least  squares  to  a fifth  order  polynomial  in  9 resulting  in  a 5 x 5 
matrix  of  coefficients  for  each  sensor  from  which  the  value  of  the 
function,  f(0,  ♦),  at  any  intermediate  value  of  9 and  ♦ may  be 
computed.  The  computer  codes  used  to  reduce  the  tabulated  function 
to  the  matrix  of  coefficients  and  to  subsequently  reconstruct  the 
function  at  arbitrary  9 and  ♦ are  given  in  Appendix  VIII.  The 
function  f ( 9 , ♦ ) itself  was  found  to  be  velocity  dependent.  This  was 
a totally  unexpected  development.  The  tunnel  was  therefore  operated 
at  its  maximum  continuous  velocity  (159  ft/sec)  and  values  of  f(9,  ♦) 
were  obtained.  The  data  were  reduced  to  polynomial  coefficients  as 
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before.  For  given  values  of  0 and  a linear  interpolation  was 
applied  so  that 


f (6,  *)  = fL(0,  <f>)  + 

V - V (59 

^ V <t>)  - f[/e»  ♦)] 

H L 

where  f^  and  f^  are  the  values  of  the  function  as  determined  from  the 
low  velocity  and  the  high  velocity  calibrations.,  respectively,  and 
and  are  the  velocities  at  which  those  calibrations  were  done. 

The  velocity  V was  determined  from  the  previous  iteration.  Use  of 
this  procedure  resulted  in  satisfactory  performance  in  test  cases. 

In  order  to  use  the  calibration  data  to  obtain  the  velocity  and 
flow  angles  from  the  three  anemometer  output  voltages,  one  first 
must  calculate  an  effective  normal  velocity  for  each  wire  by 


(E2  - E2)2 

q rs II 

K2p(T  - T )2 
s 00 

where  quantities  on  the  right  hand  side  of  the  equation  are 
for  Eq.  (56).  Equation  (58)  may  be  written  for  each  sensor 


(60) 

the  same 

as 


^1 

fl  " 

u2 

+ 

(61) 

q2 

— s 

f 2 

u? 

+ 

u3 

(62) 

q3 

at 

f3 

»? 

+ 

u2 

(63) 

fol lows . 


I 

I 

I 

I 
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In  these  equations,  is  the  value  of  the  function  f(9,  ♦)  for 
the  i1*1  sensor  which  is  evaluated  at  the  value  of  9 and  ♦ from  the 
previous  iteration.  These  equations  may  be  solved  for  the  squared 
velocity  components,  u^  . Since  there  is  no  way  to  know  whether 
Uj  is  positive  or  negative  from  these  equations,  the  proper  sign 
must  be  known  from  some  other  source.  In  the  present  experiment,  the 
orientation  of  the  probe  was  such  that  the  velocity  components  may 
always  to  assumed  to  be  positive.  Using  the  transformation  Eqs.  (47) 
through  (49),  the  velocity  components  in  the  rotor  reference  frame, 
v^ , are  computed.  The  total  velocity,  V,  and  the  resultant  flow 
angles,  0r  and  $r,  are  calculated  by  Eqs.  (64)  through  (6b). 


V = 

/vi 2 

+ v2  2 + v3  2 

(64) 

i 

0 = 
r 

sin 

(V) /V) 

(65) 

_ i 

<(>  = 
r 

sin 

(v^  /V  cos  0) 

(66) 

In  all  cases  tested,  the  method  overcorrects.  That  is,  if  the  value 
of  ♦ used  to  evaluate  the  functions  f(9,  ♦)  is  correct  but  9 is  lower 
than  the  correct  value,  9r  will  be  greater  than  the  correct  value 
of  9.  If  the  value  of  9 is  not  greatly  in  error,  $r  will  be  essen- 
tially the  same  as  the  assumed  value  of  ♦.  For  this  reason,  the 
iterations  were  uncoupled  and  under-relaxed.  Thus,  with  9 being  held 


I 

I 

I 


constant,  ♦ is  iterated  upon  using 


n 


(67) 


where  the 
♦rn  have 
the  value 


<J>n  + 1 = 0 . 8<(>n  ♦ 0 . 2 <t> 


superscripts  denote  the  iteration  level.  When  $n  and 
converged  within  the  desired  tolerance  (normally  2.5’), 
of  is  set  by  Eq.  (67)  and  0 is  updated  by  Eq.  (68), 


I 

I 

I 

I 


en+1  * 0.8  en  + o.2  e n 

r 


(68) 


When  6 has  also  converged,  the  process  is  terminated.  Note  that  when 
the  iteration  has  been  terminated,  the  correct  value  of  0 is  between 
0n  and  @rn  which  are  in  error  by  the  tolerance,  c.  Thus  choosing 
the  computed  value  of  0 to  be  the  average  of  0n  and  0fn  guarantees 
accuracy  within  c/2.  However,  due  to  the  overshoot  which  appears  to 
be  inherent  in  the  method,  allowing  the  computed  value  of  0 to  be 
defined  by  Eq.  (68)  gives  a lower  average  error.  (The  maximum 
possible  error  for  any  single  calculation  is  0.8e  when  using  this 
method  but  the  average  error  appears  to  be  on  the  order  of  0.2c.)  A 
similar  argument  may  be  made  for  the  evaluation  of  $. 
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APPENDIX  VII 


CORRECTION  OF  LEADING  EDGE  PRESSURES 


The  static  pressure  measuring  equipment  consists  of  a low  drift 
pressure  transducer  connected  to  a scanivalve,  a signal  conditioner, 
and  an  amplifier  with  an  internal  low  pass  filter.  The  average  pres- 
sures at  each  station  as  determined  by  this  system  is  added  to  the 
pressures  measured  by  the  unsteady  instrumentation  system  to  recon- 
struct the  complete  pressure  time  histories.  When  the  pressure 
coefficient  difference  for  the  leading  edge  was  plotted  (Fig.  51),  it 
became  obvious  that  something  was  wrong.  The  reduced  frequency  of 
the  distortion  is  not  high  enough  to  account  for  so  great  a shift  in 
the  data.  The  problem  was  finally  traced  to  the  amplifier  in  the 
static  data  system.  The  transducer  signal  is  first  amplified  by  the 
selected  gain  factor  and  then  filtered.  The  gain  of  the  amplifier  is 
linear  only  over  a range  of  + 10  volts  after  which  it  begins  to  clip 
the  signal.  Although  the  output  of  the  filter  never  exceeded  the 
clipping  level,  it  is  now  known  that  the  input  to  the  filter  did 
exceed  this  limit  and  was  consequently  distorted.  Thus,  the  filtered 
output  represents  the  average  value  of  the  clipped  signal.  The 
clipping  occurred  for  the  leading  edge  station  on  the  auction  surface. 
Since  this  signal  has  a relatively  large  negative  value  with  the 
variations  superimposed  on  it,  the  negative  peaks  were  clipped 


I 
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resulting  in  a slightly  more  positive  average  value  of  the  pressure 
being  indicated  than  is  actually  the  case.  The  result  is  to  shift 
the  pressure  coefficient  difference  downward  as  shown  in  Fig.  51.  No 
other  station  had  an  output  of  sufficient  magnitude  to  experience  the 
clipping  problem. 

One  expects  a pressure  coefficient  loop  at  such  a low  reduced 
frequency  to  be  centered  around  the  quasi-steady  curve.  The  loops  at 
the  other  stations  confirm  this  assumption.  Therefore,  the  loops  at 
the  two  higher  average  angles  are  plotted  with  their  centers  on  the 
quasi-steady  line.  For  the  loop  resulting  from  the  lowest  inlet 
angle,  a significant  excursion  beyond  the  static  stall  point  is 
indicated.  In  this  case,  the  portion  of  the  loop  below  the  static 
stall  point  is  centered  on  the  quasi-steady  line.  The  corrected 
loops  are  shown  in  Fig.  24. 

Now  that  the  problem  has  been  recognized,  it  is  a simple  matter 
to  monitor  the  signal  at  the  input  of  the  filter  to  detect  the  onset 
of  clipping.  If  this  occurs,  the  gain  may  be  reduced  with  only  a 
slight  loss  in  resolution.  This  monitoring  will  be  necessary  with 
testa  involving  higher  reduced  frequencies  because  the  procedure  used 
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to  correct  this  data  may  not  be  applicable. 
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APPENDIX  VIII 

COMPUTER  CODES  USED  WITH  TRIAXIAL  HOT-FILM  PROBE 

The  following  pages  are  a listing  of  a program  which  converts 
the  data  in  a table  to  a polynomial  curve  fit  in  two  variables.  A 
least-squares  technique  is  employed  to  achieve  a certain  degree  of 
smoothing  of  the  data.  The  subroutine  referred  to  as  SIMQ  is  a 
system  library  routine  for  solving  simultaneous  equations  of  the 
form 

A x = B. 

The  solution  vector  is  returned  in  B and  the  coefficient  vector 
A is  destroyed.  Any  system  routine  which  performs  the  same  function 
could  be  substituted.  In  the  event  that  one  is  not  available,  a 


listing  of  SIMQ  is  included. 
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PROORAM  FOR  GENERATING  HAST  SQUARES  FIT  TO  DIKEC1IONAL  CALIBRATION 
DIMENSION  F*3.9.9).THETA<20>.PMI  ( 20  > . A<  3 . 9 . 5 1 . YDUM*  20 ) .COEF  ( 20  i 
RTD-37. 293/8 

READ  IN  F (NWIRE . ITHETA.JPMI ) AND  THE  THEIA  AND  PHI  ANN AT b 

READ* 1 .900) << THETA* I > .PHI < J ) . < F * N. I . J ) .N- 1 . 3 > • J* I .9 > • I * 1 .9 ) 

DO  100  I JR- 1.9 

THETA* I JR  > -THETA* 1 JN  > /RTD 

PHI  < 1 JR  >-PHI  ( I JR) /RTD 

THE  FOl  LOWING  SUBSCRIPT  NOTATION  WILL  BE  USED 
N HIRE  NUMBFR 
I THETA  SUBSCRIPT 
J PHI  SUBSCRIPT 
R I R-TH  HARHUN I C IN  THEIA 
RJ  R-TH  HARMONIC  IN  PHI 
ALL  OTHERS  ARE  DUMMY  SUBS*  KU  IS 

FIT  F AS  A FUNCTION  Of  PHI  FOR  EACH  WIRE  (N)  AND  THETA  (I) 

DO  300  N-l  . 1 
DO  300  I - I • 9 

PUT  F INTO  A DUMMY  VFI  FOR 
DO  200  J- I ,9 
YDUM  < J)  -F  <N.  I .J) 

CONFUTE  POLYNOMIAL  I I AS  I SOUARt  S MI 
CAL  L POl  YS(I<  PHI  . YDUM  •CUE  f . 9 RE  I > 

STORE  THESE  (OFFFICIENTS 
DO  300  RJ-1.5 
A < N»  I • R J > - (.  UE  F <RJ> 

CURVE  FIT  THESE  COEFFICIENTS  AS  FUNCTIONS  OF  THETA 
DO  300  N- 1 » 3 
DO  300  R 1-1.3 
DO  400  1-1.9 
YDUM* I >-A<N.I ,RJ) 

CALL  POLYSO* THETA. YDUM. COEF .V.3.RRET) 

STORE  FIT 

DO  300  R 1-1.3 
A < N • R I .RJ) "COEF ( R I ) 

COEFFICIENTS  ARE  NOW  CALCULATED...  WRITE  INTO  DATASET. 

DO  600  N-l.  3 

WRITE *4. 901 ) < * A*N. I . J) . J- 1 ,3) .N. I. 1*1 .5) 

STOP 

F ORHAT  < 3F 1 3 . 7 ) 

F0RMAT*3E14. 7.215) 

END 

SUBROUTINE  POLYSO*  X. Y.B.NPTS.NTERMS.KRFT ) 

POLYNOMIAL  LEAST  SOUARFS  FIT 

X INDEPENDENT  VARIABLE 

Y DEPENDENT  VARIABLE 

B COEFFICIENT  VECTOR  FOR  POLYNOMIAL 

NPTS  NUMBER  OF  POINTS  TO  BE  FIT 

NTERMS  NUMBER  OF  TERMS  TO  BE  IN  POLYNOMIAL 

KRET  0 INDICATES  GOOD  FIT...1  INDICATES  FIT  FAILED 

DIMENSION  X<20) .Y(20) .B<20) .A* 20.20) 

CHECR  FOR  AN  APPROPRIATE  ORDER  OF  EQUATION 
IF < NTERMS. GE. NPTS > NTERMS^ NPTS- 1 
FORM  MAR T I X OF  COEFFICIENTS 
DO  200  I -1. NTERMS 
DO  200  J-l. NTERMS 
IF ( IF J.NE *2>G0  TO  100 
A*  I .J) -FLOAT* NPTS) 

00  TO  200 

A* I. J)-0.0 

IEXP-I+J-2 

DO  199  K- 1 . NPTS 

A*  I. J)-A(I. J)FX(K)**IEXP 

CONTINUE 


00000010 
00000020 
'U/OOOO  JO 
00000(/40 
00000030 
00000060 
0000007 o 
OOOOOOBO 
00000090 
00000100 
00000110 
00000 1 20 
OOOOOl 20 
00000140 
OOOuv.  1 30 
00000160 
OOOOOl 70 
OOOOOl bo 
OOOOOl VO 
000o0200 
000002 lo 
00000220 
00000230 
00000240 
00000230 
00000260 
00000270 
00000280 
00000290 
00000300 
OOOOOl 1 0 
00000320 
00000330 
00000340 
00000350 
00000360 
00000370 
00000380 
00000390 
00000400 
00000410 
00000420 
00000430 
00000440 
00000450 
00000460 
00000470 
00000480 
00000490 
00000500 
00000510 
00000520 
00000530 
00000340 
00000550 
00000560 
00000570 
00000580 
00000590 
00000600 
00000610 
00000620 
00000630 
00000640 
00000630 
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C FORH  VECTOR  FOR  RIGHT  HAND  SIDE  OF  EOUATION  00000660 

•< 11-0.0  00000670 

DO  300  K-1.NPT8  00000680 

300  D< 1 >-D< 1 )4Y<K)  00000690 

DO  400  I-2.NTERHS  00000700 

D<I>-0.0  00000710 

IEXP-I-1  00000720 

DO  400  K-1.NPT9  00000730 

400  D(I)«D(I>4Y(K>«X(K>*tIEXP  00000740 

CRLL  8IMKR. D. 20. NTERHS.KRCT ) 00000750 

H (KM  T .((i.MM  TURN  OOOvv 

VKITt  (3.500)  0..0«07  70 

500  HlKMATl/  MIIMMiMtlt*  SIMUUIAR  SOI  111  ION  ,ou  on 

Kf  TURN  Ouvlru'-V'j 

^ HO  OCOOobuO 

END  Of  DATA 
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The  polynomial  coefficients  may  be  used  to  reconstruct  the 
function,  f,  at  any  arbitrary  value  of  theta  and  phi  by  the  algo- 
rithm listed  below.  Since  this  algorithm  will  return  a result  for 
values  of  the  angles  far  outside  the  range  of  the  original  data 
(and  this  result  may  well  be  in  error  because  the  polynomial  fit  is 
not  a good  means  of  extrapolation),  one  must  be  careful  in  applying 

FUNCTION  FPT  (PHI,  THETA,  N) 

COMMON/ DCOF/A(  3,5,5) 

FPT-0 . 0 
DO  200  1-1 ,5 
C“A(N,5 ,6-1) 

DO  100  J-2,5 
C-C*THETA+ A(N,6-J ,6-1) 

FPT«FPT*PHI+C 
RETURN 
END 

i ■ § |i]j| 

t 

i; 

I! 


100 

200 
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SUBROUTINE  SI MO  < A » B » ND » N » KS ) 

DIMENSION  A ( ND » ND > fB(ND) 

LOGICAL  *1  BIGEXF' » DI VEXF' 

EQUIVALENCE  ( DI  V r D I VEXF' ) » < BIGA » BIGEXP ) 
REAL  DIV/ZOOIOOOOO/  » T0L./Z3C800000/ 

DO  300  1 = 1 »N 
BIGA--0. 

DO  100  J= 1 » N 

100  BIGA--AMAX 1 ( BIGA f ABS ( A ( I f J ) ) ) 

DI  VEXF'=BIGEXF' 

DO  200  J=1  f N 
200  A( I f J>=A( I f J)/DIV 
300  B(I)-B<I)/DIV 
KS  = 0 
JJ=-N 

DO  9 J= 1 f N 
JF'l-- J+l 
B I G A = 0 . 0 
DO  2 I = J r N 

IF  < ABS  (BIGA)- ABS  ( A(  I f J)  ) •)  1 f 2 » 2 

1 BIGA=A ( I f J) 

I MAX  = I 

2 CONTINUE 

IF ( ABS ( B IGA ) -TOL ) 3 r 3 f 4 

3 KS  = 1 

IF (BIGA.EQ.O. ) RETURN 

4 DO  6 K = JrN 
IF< IMAX-J>5f6f5 

5 SAVE=A  < I MAX » K ) 

A ( I MAX » N ) —A ( J f K ) 

A ( J f K ) = SAVE 

6 A(  JfK>=A<  JrlO/BIGA 
SAVE-  Ft  ( I MAX ) 

B ( IMAX ) =B ( J ) 

B(J)=SAVE/BIGA 
IF< J-N)7f 10 f 7 

7 DO  9 IX=  JF'l  r N 
DO  8 JX-JF'  1 f N 

8 A ( I X f JX  ) =A  ( I X f JX ) -A  ( I X f J ) * A ( J r JX  > 

9 B< IX)=B( IX)~B( J>*A< IXf J) 

10  NY  = N-  1 

DO  11  J- 1 r NY 
IB- N- J 
IC  = N 

DO  11  K=1fJ 

B(IB)-B(IB)-A(IBfIC)#B<IC) 

11  IC=IC-1  '■ 

RETURN 

END 


I 
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